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Abstract. Active galactic nuclei (AGN) at the center of galaxy clusters with gas cooling times that are much shorter than the 
Hubble time have emerged as heating agents powerful enough to prevent further cooling of the intracluster medium (ICM). 
We carried out an intensive study of the AGN heating-ICM cooling network by comparing various cluster parameters to the 
integrated radio luminosity of the central AGN, L R , defined as the total synchrotron power between 10 MHz and 15 GHz. This 
study is based on the HIFLUGCS sample comprising the 64 X-ray brightest galaxy clusters. We adopted the central cooling 
time, f coo i, as the diagnostic to ascertain cooling properties of the HIFLUGCS sample and classify clusters with / coo i < 1 Gyr 
as strong cool-core (SCC) clusters, with 1 Gyr < f cocJ < 7.7 Gyr as weak cool -core (^VCC) clusters and with t coo \ > 7.7 Gyr as 
non-cool-core (NCC) clusters. We find 48 out of 64 clusters (75%) contain cluster center radio sources (CCRS) cospatial with 
or within 50 kpc of the X-ray peak emission. Furthermore, we find that the probability of finding a CCRS increases from 
45% to 67% to 100% for NCC, WCC, and SCC clusters, respectively. 

We use a total of ~ 140 independent radio flux-density measurements, with data at more than two frequencies for more 
than 54% of the sources extending below 500 MHz, enabling the determination of accurate estimates of Lr. We find that Lr 
in SCC clusters depends strongly on the cluster scale such that more massive clusters harbor more powerful radio AGN. The 
same trend is observed between Lr and the classical mass deposition rate, Mciassdcai in SCC and partly also in WCC clusters, 
and can be quantified as Lr oc M'^f 1 ^ 25 . We also perform correlations of the luminosity for the brightest cluster galaxy, Lbcg, 
close to the X-ray peak in all 64 clusters with L R and cluster parameters, such as the virial mass, M 50 o, and the bolometric 
X-ray luminosity, Lx- To this end, we use the 2MASS /f-band magnitudes and invoke the near-infrared bulge luminosity-black 
hole mass relation to convert Lbcg to supermassive black hole mass, Mbh- We find a weak correlation between Mbh and Lr for 
SCC clusters, L R ~ M^ >±0A2 , although with a few outliers. We find an excellent correlation of L B cg with M 50Q and L x for the 
entire sample, the SCC clusters showing a tighter trend in both the cases. We discuss the plausible reasons behind these scaling 
relations in the context of cooling flows and AGN feedback. 

Our results strongly suggest an AGN-feedback machinery in SCC clusters, which regulates the cooling in the central 
regions. Since the dispersion in these correlations, such as that between Lr and M c iassicai or Lr an d Mbh, increases in going from 
SCC to WCC clusters, we conclude there must be secondary processes that work either in conjunction with the AGN heating 
or independently to counteract the radiative losses in WCC clusters. 



1. Introduction 

In recent years, heating by active galactic nuclei (AGN) 
through outflows has gained fundamental importance in the 
realm of la rge-scale structure and galaxy formation. Several 
studies (e.g.lSiiacki et al.|2007llCroton et al.l2006l:lBower et al 



2006t IScannapieco et alll2005t ISilk & Reeslll998l) have con 
tributed to this comprehensive picture, wherein AGN feedback 
is considered an attractive solution to several connected prob- 
lems, suc h as the high-mass end truncation of galaxy distribu- 



tion (e.g. lBenson et al. 2003) and the absence of cooling -flows 



ing- 

:n2 



in centers of galaxy clusters (e.g. IMcNamara & Nulsenll2007l) . 
According to these studies, AGN heating at the centers of clus- 
ters may likely be responsible for quenching condensation of 



the hot intracluster medium onto the cluster galaxies, thereby 
giving rise to the cutoff at the bright-end of the galaxy luminos- 
ity function and also regulating the cooling flows. 

Gas in the ICM cools via X-ray emission. In the centers of 
some clusters, the high density leads to significant loss of en- 
ergy, such that the gas radiates away all its energy in a short 
(<k 1 ///o) time. In the absence of any heating mechanisms, 
in order to support the overlying gas and restore hydrostatic 
equilibrium, there is a steady inflow of gas towards the cluster 
center, which is ofte n referred to as the classical cooling flow 
model (lFabianlll994l) . These so-called cool-core clusters (CC) 
have centrally peaked X-ray surface-brightness profiles imply- 
ing gas cooling times orders of magnitude shorter than the age 
of the cluster. However, (1) the high resolution XMM-Newton 
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RGS spectra of CC clusters have not found the expected 
amounts of cool ga s in their cores (e.g . Tamura et al 



2003 



Peterson et"ai]|200 it iKaastra et alj|200lfc [Peterson et al. 
Xu et al.ll200alSakeluou et al.l 120021: 1 Sanders et al.ll2008l) . and 



2001 



(2) even though the cooling of the ICM is manifested in 
the form of on-going star formation obse rved in the bright 
est cluster galaxy of several clusters (e.g. iMittaz et al.l 12001 



Alien! 1 19951) . it is far below t he predicted amount of the star 
formation rates and CO (e.g. iMcNamara & O 'Connelll[l989l 
Edge & Frayerll2003l) . Additionally, the gas temperature in the 
central regions as determined from the X-ray spectra of these 
clusters is much higher than that expected based on the cooling 
flow model and has been fou nd to drop not muc h below 40% 
the ambient temperature (e.g. lHudson et al ] |2008h . 

Several heating strategies have been proposed to overcome 
the cooling flow problem. Feedback from supernovae is an 
important form of heating but has been shown to be suffi- 
cient to balance energy losses only in lo w-luminosity ellipticals 
with shallower gravitational potentials dMathews & Brighenti 
20031) . 

Another heating scenario is thermal conduction which 
leads to an inward heat flow fr om the outskirts of the galaxy 
clusters. IVoigt & Fabianl (120041) have shown that even though 
thermal conduction may provide enough heating to offset cool- 
ing in the hotter (T>5 keV) part of the clusters, the central parts 
of the cooling region remain largely unaffected by this pro- 
cess. Similar to supernovae heating, thermal conduction also 
has the effect of only slowing down the evolution of intraclus- 
ter mediu m by causing the cooling time to increase by a factor 



2005) but leaves the cooling catastrophe 



of a few dPope et al 
inevitable. 

In this work, we focus on the self -regulated AGN feedback 
as the current favored mechanism to expla in the dearth of 
coolin g by-products in galaxy clust er s (e.g . Voit & Donahue 

20021: 



20051: Rovchowdhurv et al. 2004 : Churazov et al 



Binnev & Taborl 1995 ). In this framework, accretion of the 



cool collapsed intracluster medium (ICM) ignites the central 
active galactic nucleus, which returns a fraction of the accreted 
power back to the ICM. The bulk of the energy transfer is 
believed to happen through mechanical dissipation of the AGN 
power. The lead evidence comes from the observations of 
numerous galaxy clusters featuring X-ray deficit low density 
regions, known as cavities. Such cavities have been observed 
to correlate spectacularly with radio jets and lobes indicating 
that they are likely regions emptied of ICM by the e xpanding 
rad io lobes [e.g. Pe r seus. iBoehringer et alj (19931): Hydra - 
A, IMcNamara etal] d200oi): A 2052, iBlanton et all d200lh : 



A2597 IMcNamara et al.l d200ll); A4059. iHeinz et all J2OO2I) : 
ySui 



A478. ISun et all (120031) : A2029. IClarke et al.l d2004l) : A2199 



Gentile et al.l d2007l) l. The AGN-blown cavities transfer heat 



to the I CM potentially by generating sound and weak shock 



2002; Fabian etal. 2003; Mathews et al 



wave s dJones et al 
20061) . by doing pdV work against the ambient medium and 



dissipation of cavity enthalpy in t he wa k e of buoyantly risin g 



cavities (e.g 
In addition to 



Ruszkowski et al. 2004; Birzan et al. 2004) 



direct AGN mechanical heating via radio 
bubbles, there are also alternative proposed mechanisms such 
as AGN cosmic-ray heating combined with convection (e.g. 



Chandran & Raseral 12007k or conduction (e.g. iGuo & Ohl 
2008tlVoitetalJl2008l) . 



Numerous results over the last couple of decades have 
confirmed that radio-loud AGN dwell preferentially in bright- 
est group and cluster galaxies (B CGs), as opposed to other 
galaxies of the same stellar mass ( von der Linden et all 2007 



Best et alJl2007T:lBagcM & Kapahilll994tlValentiin & Biileveld 



19831) . It has also been found in these and other studies that the 
CC clusters are particularly conducive for cD galaxies which 
are radio-loud, even though the fraction of radio-loud cD galax- 
ies in CC cl usters varies from study to study ranging from 
70% to 95% dBurnsI 1 990tlEdwards et al J2007t [Dunn & Fabian 



2006). The spread in the fraction between different studies can 



be attributed to the varying selection criteria used for con- 
structing cluster samples and the use of not so up-to-date X- 
ray and radio observations. The latter effect may result in the 
same cluster being identified as a CC cluster in some works 
and a non-cool-core (NCC) cluster in others. A few examples 
being A 1650, which based on Einstein o bservations has been 
quoted as a NCC cluster by Burns ( 1990l) . but w hich our data. 



based on high-resolution Chandra observations dHudson et al 



120081) . clearly reveal it to have a cool core with a central cool- 
ing time of about a gigayear and a pre dicted mass depositio n 
rate of about 100 M Q yr _1 (also see iDonahue et al.l 120051) . 



Similarly, A3 158 and A3 195 have been identified based on 
lo w-sensitivity and low -resolution ASCA data as CC clusters 
bv lEdwards et al.l (12007b . Our results imply otherwise; both are 
merging systems each with a central cooling time longer than 
12 Gyr and the expected mass deposition rates being consis- 
tent with zero. Despite these inconsistencies, most studies are 
by and large in agreement with one another and set the average 
abundance of radio-loud CC clusters in the local Universe to 
around 80%. 

Recent analyses of galaxy clusters have shown that of those 
CC clusters which require heating, at least 40% harbor cavities 
that contain sufficient energy to balance the radiative losses 
dRaffertv et all l2008t iNulsen et al.lE006l iRafferty etaD 12006 : 
Dunn & Fabianl2006 ). However, the details of the various heat- 



ing mechanisms set into motion by the central AGN are not 
clear and are issues currently under investigation. Also, the 
local conditions in the ICM that lead to a quasi-steady state 
of gas deposition onto the central regions, and presumably 
onto the supermassive black hole, and the concomitant AGN 
heating of the ICM either periodically or continuously remain 
largely unknown. A parallel model that is emerging to ex- 
plain the deviation of the observed cluster properties, especially 
the entropy profiles, from the predictions of the pure cooling 
model is linked to preheating or entropy injection at incipi- 
ent stages of cluster fo rmation, even prior to cluster collapse 
dMcCarthv et al.l 120071 and references therein). Preheating is 
entailed by reduction of central densities, hence, central lumi- 
nosities, which leads to flat cluster entropy profiles. This mod- 
ification along with post cluster-formation processes, namely 
radiative cooling and gravitational heating, provides a better 
match to the observed entropy profiles of galaxy clusters. Yet, 
while preheating may alone account for the differences in en- 
tropy profiles in NCC clusters, catastrophic cooling at small 
radii in CC clusters can still not be bypassed. In order to main- 
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tain their entropy profiles at observed levels, one or more addi- 
tional sources of on-going heating are required. 

In this work, we aim to gain a more comprehensive under- 
standing of the AGN-regulated cooling and heating. We scruti- 
nize the ways in which AGN heating is connected to the cooling 
of the ICM based solely on the total radio (synchrotron) output 
of the AGN. To achieve this goal, we use a sample of galaxy 
clusters for which there exist complete radio and X-ray data. 
We take our analysis further by examining the scaling relations 
between the BCG near-infrared luminosity and cluster param- 
eters (mass and luminosity). We derive the mass of the super- 
massive black hole using the near-infrared bulge luminosity- 
black hole mass relation and inspect whether there is a relation 
between the black hole mass and AGN radio luminosity. These 
correlations are made taking into account the possibility of ob- 
taining different relations depending on the cool or non-cool 
type cluster environment. The improvement over previous anal- 
yses lies in the quality of the cluster sample and of the available 
X-ray and radio observations. 

We describe the sample in Section |2j giving details about 
radio and X-ray data and related quantities in Section |2~T1 and 
Section l2~2l respectively. The results are presented in Section[3] 
including fractions of CC and NCC clusters with and without 
central radio sources in Section ITT] cooling and AGN activity 
in Section [3T2l and correlations of the BCG luminosity with ra- 
dio and X-ray parameters in Section [3~3l We discuss our results 
in Section|4]and end with conclusions in Section[5] Throughout 
this paper, we assume the ACDM concordance Universe with 
H = A 7 i71 km s" 1 Mpc 1 , O m = 0.27 and Q A = 0.73. 

2. Our sample 

The goal of this study is to cross-correlate the cooling activity 
with the presence of a radio galaxjQ cospatial with the peak 
of the cooling flow region or, synonymously, the peak of the 
X-ray emission. Further, we want to investigate whether there 
is a special coupling between the AGN activity and its sur- 
rounding hot cluster medium in CC clusters as compared to 
NCC clusters. This calls for an objectively selected sample. 
Samples selected based on what is available in public archives 
are subject to unknown selection effects ("archive bias"). For 
example, the fraction of CC clusters in the Chandra archive 
may be biased higher (or lower) than the same fraction of 
clusters in the XMM-Newton archive. To this end, we con- 
duct our study b ased on the largest X-ray flu x-limited sample, 
the HIFLUGCS dReiprich & Bohringerll2002l) sample, selected 
from the ROSAT All-Sky Survey outside the Galactic plane us- 
ing the flux limit, f x (0.1 - 2.4) keV> 2 x 10~ n erg s" 1 cm" 2 . 
This sample comprises the 64 X-ray brightest clusters and 
spans a redshift range 0.0037 < z < 0.2153 with the mean 
(z) ~ 0.05. All 64 clusters have been observed with Chandra 
and all but one (A2244) have been observed with XMM-Newton 
to acquire high quality X-ray data. In this paper we make use 
only of the Chandra analysis because we are interested in the 
cluster cores and Chandra has currently the best spatial reso- 



lution (0. 5) of any X-ray telescope. We note that since flux- 
limited samples are as such biased towards clusters with high 
or boosted luminosities, this sample might seem to favor CC 
clusters preferentially over NCC clusters. At any given redshift, 
CC clusters are more likely to be picked up than the NCC clus- 
ters due to their enhanced central luminosities. However, merg- 
ing clusters present the same bias as cluster merger events en- 
tail tempo ral enhancements in the global luminosities and tem- 
peratures (iRicker & Sarazinl 1200 ll). Since merging clusters ar e 
mostly NCC clusters fo'Hara etd]|2006tlHudson et alJEooi) . 
this may balance out the former selection bias towards, at least, 
the strong (to be defined in Section 13.11 ) CC clusters. We 
note that even in the presence of a bias against transition clus- 
ters (neither strong CC clusters nor NCC), objectively selected 
samples, such as HIFLUGCS, can be directly compared to sim- 
ulated flux-limited samples, at both low-z and high-z, and the 
bias may be calculated. 

2.1. Radio data 



We compiled and in many cases reanalyzed radio observations 
of all 64 clusters from either literature or archives (VIA, NVSS, 
VLSS and MOST) to study the radio properties of the centrally 
located AGN in the HIFLUGCS clusters. The data from the 
archives were processed in the Astronomical Image Processing 
Software (AIPS) package provided by NRAO. 

High-frequency archival radio data (500 MHz) were an- 
alyzed using the standard data-reduction procedures within 
AIPS, wherein the resulting map usually constituted of a single 
hybrid image. Low-frequency data (330 MHz and 74 MHz), in 
addition, were carefully analyzed to remove bad data affected 
by radio frequency interference (RFI) using the AIPS tasks, 
SPFLG and TVFLG. In case of pseudo-continuum mode ob- 
servations, the effects of bandwidth smearing were tackled by 
keeping the data separate over the spectral channels. And lastly, 
in order to correct for 3D effects and image degradation due to 
bright sources far away from the phase-center, we employed 
the 3D-imaging feature embedded in the AIPS task, IMAGR. 
This technique of 3D imaging results in a mosaic of "facets", 
each of which is independently and simultaneously "cleaned" 
to yield a final map of the entire primary beam. 

2.1 .1 . Location of a cluster radio source 

Before pursuing the cross-correlation of the global X-ray prop- 
erties of our sample with the radio properties of cluster central 
radio sources (CCRSs), we need to establish a criterion based 
on which to identify central radio sources. Several works have 
presented evidence for a special relationship between the cool- 
ing activity in cluster cores and BCGs located within a cer- 
tain distance to the X-ray peak. Recently. Edwards et al.l (12007b 
found in their study based on two samples, the NFPS data set, 
an X-ray selected sample and the C4 catalog, an optically se- 
lected sample built from the SDSS, that only those BCGs that 
lie within 70 hz\ kpc of the X-ray peak of a cooling flow clus- 



1 We use the terms "radio galaxy" and "AGN" interchangeably 
throughout this paper. 



ter have significant line emission. Even though the optical line 
emission observed in these BCGs can be inferred either as 
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Fig. 1. The separation between the BCG and the X-ray Peak. Plotted 
on the X-axis are numbers assigned to each cluster arbitrarily. The 
filled circles (blue) represent BCGs which harbor a radio source and 
the crosses (red) correspond to BCGs without one. The black horizon- 
tal line at y = 50 kpc corresponds to the criterion for determining 
whether or not a cluster has a CRS and the grey horizontal line at 
y = 12 Kj\ kpc corresponds to the worst uncertainty associated with 
the position of the X-ray peak in the sample. 



a sign ature of AGN activity or star formation, Edwards et al 
(12007b also show that 74% of the strongly emitting BCGs in the 
SDSS sample, defined as having the H a equivalent width > 2A, 
have the diagnostic emitting-line ratios characteristic of AGN 
activity (with a likely higher fraction for the BCG). Therefore, 
the formal basis of marking a radio source as "central" in the 
study presented here was to have the AGN within 50 hz\ kpc 
of the X-ray peak, a mor e conservative limit than that proposed 
by Edwards et alJ d2007h . It is noteworthy that this cut (as op- 
posed to a more stringent cut of 12 kpc, see below) had 
actually to be invoked only for four clusters. These four clus- 
ters are A3562, A2142, A4038 and A3376 with the X-ray peak 
and BCG separation as 30.4 h~\ kpc, 21.8 Kj\ kpc, 14.9 h~\ kpc 
and 14.2 kpc, respectively. For the rest of the sample, the 
flagging was straight forward in that the separation between the 
X-ray peak and the radio active BCG was less than 12 hz\ kpc. 
The 12 hz\ kpc as the yardstick comes from the fact that since 
the HIFLUGCS clusters span two orders of magnitude in red- 
shift, the Chandra resolution implies varying accuracies with 
which the X-ray peak may be determined for different clusters 
and 12 kpc corresponds to the worst 1-cr uncertainty. The 
separation between the BCG and the X-ray peak for all 64 clus- 
ters is shown is FigureQ] Applying this criterion, we find a total 
of 48 clusters with centrally located radio sources. 



2.1 .2. Integrated radio luminosity of a CCRS 

One of the primary concerns in accumulating the radio data 
was to have a good spectral coverage, particularly, at the low- 
end of the radio spectrum. Of the 48 CCRSs, 65% have radio 
data below 500 MHz and 46% have radio data below 80 MHz. 
Low-frequency radio observations are important to perform a 
full spectral analysis for these clusters for two reasons. First, 



since the main contribution to the synchrotron radio luminosity 
comes from the low-end of the radio frequency spectrum, this 
will result in precise determination of the synchrotron power 
in these systems. In our subsequent study, the energy in the 
radio emitting particles will be compared to the mechanical 
energy of the cavities, and thereby the partitioning of energy 
between radiation and mechanical (cavity) energy can be mea- 
sured. Second, with dense enough spectral sampling, spectral 
breaks may also be visible (as have already been seen in a 
few CCRSs in the HIFLUGCS sample, see below). A spectral 
break in a system with cavities is an extremely useful observ- 
able since it is representative of the time since the last injection 
event or particle production, and therefore, a good indicator 
of the age of the cavity emission. The CCRSs, which we al- 
ready know to show spectral breaks are presently being cross- 
checked with presence of cavities for a future study. 

For sources with no observations or confirmed detection be- 
low 100 MHz, the low-end of their spectra were constrained us- 
ing the 74 MHz VIA Low-frequency Sky Survey (VLSS) with 
an average point-source detection limit of 0.5 Jy/beam and a 
resolution of 80 ". Shown in the left panel of Figure [2] are the 
spectra of all CCRSs. Shown in the right panel of Figure [2] 
are the spectra of only a subset of CCRSs that show interest- 
ing features such as spectral breaks and turn-overs indicative 
of spectral aging and synchrotron self-absorption, respectively. 
Amongst the CCRSs shown in the right panel, we note there 
are also a few clusters which show spectral steepening at low- 
frequencies, which we believe is due to the superposition of dif- 
ferent radio components with different spectral properties (due 
to varying sizes and distances from the central engine). 

The details of the radio data used for this work, such as the 
various frequencies and the corresponding flux-densities used 
to estimate the radio luminosities of the cluster central radio 
sources, along with the references can be found in TableQ] The 
synchrotron radiation is assumed to have a powerlaw spectrum 
given by S (v) oc v~ a , where S (v) is the flux density at frequency 
v. Thus, the integrated rest-frame radio luminosities, Lr, of the 
CCRSs were calculated by step-wise integration: 



L, + i = AnD: 



]So I (v 



4-nDj- 
' (1 



S v 



■ <Xi+\,i) 



(I) 



where So is the flux density of the radio source at either of 
the two rest-frame frequencies, v, or v, + i, L,+ i is the radio lu- 
minosity in the frequency range [v,-, v,+i], aj + u is the spectral 
index between v, and v,+i, and D\ is the luminosity distance. 
The total radio luminosity was calculated by extrapolating the 
spectral indices obtained at the lowest observed frequency to 
10 MHz and at the highest observed frequency to 15 GHz. 
Thus, L tot = 2 Li+\- Of all the clusters with CCRSs, 27 have 
(reliable) data at more than two frequencies and 18 have data 
at two frequencies. The remaining three clusters (A576, A3 158 
and A3562) have data at only one frequency and we used a — 1 , 
the average spectral index of the CCRSs in our sample, to cal- 
culate their total radio luminosity. 
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Frequency (GHz) Frequency (GHz) 

Fig. 2. The radio spectra of the cluster central radio sources. Shown in the left panel are the spectra of all 48 CCRSs, where the symbols are 
the actual measurements. Most of the measurements at 74 MHz and 1.4 GHz are taken from VLSS and NVSS, respectively. Shown in the right 
panel are the spectra of only those CCRSs which show spectral curvature. The black horizontal solid line represents the VLSS point-source 
sensitivity constraint at the VLSS observing frequency of 74 MHz indicated by the black vertical line. 



Even though the formal errorbars for Lr were derived us- 
ing the background root-mean-square in the maps, these do 
not take into account the uncertainty arising due to the lack of 
knowledge of the shape of the radio spectra down to the lowest 
frequencies, except for a handful of radio sources well-studied 
at all radio-frequency bands (such as, Hydra-A (A0780), 
Centaurus (A3526), A1795, A2029, A2052, A2199, A2597 
and A4059). Radio sources often show a spectral turn-over 
at low-frequencies. This is attributed to synchrotron self- 
absorption which kicks in with increasing optical depth and is 
manifested by a rising spectrum with a = -2.5. Two examples 
of spectral turn-overs can be easily seen in the right panel of 
Figure|2] Considering the possibility of other CCRSs showing 
similar turn-overs, we calculated a lower-limit on Lr based on 
the assumption that the spectra of the CCRSs turn over right 
below the lowest observed frequency. The difference between 
the lower limit derived in this manner and the integrated ra- 
dio luminosities assuming the continual of the spectra beyond 
the lowest observed frequency provides a more realistic, albeit 
conservative, uncertainty on Lr. This is the reason for having 
highly asymmetric errorbars for Lr, as can also be seen in many 
of the plots. 

In other works (e.g. IBurns Il99d: iPeres et al.lll998l) . very 
often the comparisons between the quantities representing the 
cooling flow strength in clusters and the radio power of the 
BCGs are based solely on the monochromatic radio luminosi- 
ties, such as, the 1.4 GHz luminosities derived from NVSS or 
FIRST or 5 GHz luminosities derived from the Green Bank 
Survey. A useful exercise followed up with our data was to 
compare the integrated radio luminosities to L l 4 GHz> where 
the latter is given by 



l ia GHz = 4nD J s 1.4 GHz(! + Z> 



ati.4,,—1 



(2) 



where a\Aj is the spectral index between 1.4 GHz and 
one of the two neighboring frequencies. Shown in Figure [3] is 
the integrated radio luminosity versus the monochromatic ra- 
dio luminosity at 1.4 GHz. Also shown is the best-fit power- 



law derived using the FITEXY least-squares line-fitting routine 
( Press et aDl 1 992 ) given by 



lO 42 /^ 2 ergss- 1 



= a x 



A A GHz 



10 32 /i7 2 ergss- 1 Hz" 1 , 



(3) 



where a = 1.04 + 0.03 and b = 0.98 + 0.01. This algorithm al- 
lows for fitting in only one parameter (Lr in this case) but takes 
uncertainty in both X and Y into account. Our study shows that 
there is a fairly good correlation between the Lr and L x 4 qhz- 
This is not surprising since once a pedestal value for radio 
sources is determined, the total power should scale with the 
spectral index, which for the CCRSs in our sample is quite sim- 
ilar and centers around unity. There are, however, a few CCRSs 
(such as 2A0335, A3376, MKW3S and A4038) which are in- 
consistent with the best-fit relation, all of which have spectral- 
indices steeper than unity at low frequencies and were excluded 
while determining the best-fitting powerlaw. Hence, we con- 
clude that even though the monochromatic luminosity (in this 
case at 1.4 GHz) is a good proxy for the total radio luminosity 
and may be used in cases where additional spectral information 
is not available, for precise radio correlations demanding least 
intrinsic scatter, the total radio power computed from detailed 
spectral analyses should be used. 



2.2. X-ray data 

The complete HIFLUGCS sample has Chandra observations 
and the data have been homogeneously reprocessed using 
CIAO 3.2.2 and CALDB 3.0. For a detailed description of the 
X-ray data -reduction, and data- and error-analysis, the reader is 
referred to lHudson et al.l (120081) . In the following we briefly de- 
scribe only those cluster X-ray parameters that are meaningful 
in the context of the present work. 
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Fig. 3. A comparison between the radio luminosity integrated be- 
tween 10 MHz and 15 GHz and the monochromatic radio luminosity 
at 1.4 GHz (red solid line). The dotted grey lines correspond to l-<x 
deviation above and below the best-fit line. 




Fig. 4. The central cooling-time distribution clearly showing the steep 
rise in the fraction of clusters with cooling times shorter than 1 Gyr. 
For clarity, we show the distribution in both linear as well as log scale 
(inset panel) on the x-axis. 



2.2.1. Cluster parameters 

We extracted the central density profiles by fitting either a sin- 
gle or double [}- model to the surface-brightness profile annuli. 
Similarly, the central temperature profiles were derived by fit- 
ting spectra to annuli to an absorbed thermal model. From the 
temperature and density profiles, the central cooling time at 
r = 0.0047?5oo (0.4% /?5oo), fcooi, entropy, K, and cuspiness, 
a, were calculated as per the following expressions: 



^cool 



5 (WjQ + »eo)fcr 48 

2 « e 2 A(r 48 ) 



K(r) = kT(r)nf /3 (r) 

d\og(n e ) 

a — 

dlog(r) 



where «io and « e o are the central ion and electron densities, 
respectively, determined at r — 0.004T?5oo, 7 48 is the average 
temperature of the - 0.048i?soo region, A{T^) is the cooling 
function for a plasma at T/& and n e (r) is the electron density at 
a given radius r from the cluster center. The cuspiness is calcu- 
lated at a distance r = 0.04J?5oo from the cluster center. Here, 
7?5oo is the radius within which the average cluster mass density 
is 500 times higher than the critical density of the Universe. 

To derive the cluster entropy profiles, the best fit density 
profiles were binned in steps of 2 . Then for each bin, the 
value of temperature corresponding to that radius was adopted 
to calculate the entropy for that bin using the expression in 
Equation |U Since the annuli created to derive the temperature 
profiles need not necessarily coincide with the 2" density bins, 
in the case where there was a jump in the temperature within 
a density bin, the average value of the two temperatures was 
used. 

The virial temperature of the cluster, used as a scaling pa- 
rameter in some of the cross-correlations presented in the forth- 
coming sections, was determined by fitting the temperature 
profiles to broken powerlaws. This was done so as to prevent 
the cool-core gas from biasing the estimate for the global 'viral' 



temperature, 7\,;, .. In those cases, where the inner powerlaw had 
a declining slope towards the center, which is representative 
of the cool gas at the centers of CC clusters, we excluded the 
core region as determined from the break in the powerlaw from 
the fit. From the estimate of T wn , the virial mass, M500, within 
i? 500 w as determined using the relation by iFinoguenov et al 



feOOlh . M500 = akT^ 1O I3 M , where a = 2.5 + 0.2 h^ 1 and 
b = 1.676 ±0.054. 

The classical mass deposition rates, Mciassicai, were derived 
from the gas temperature and density profiles. M c i ass ic a i is the 
ratio of the total gas mass within the cooling radius, r coo i, 
defined as the region at which the ? coo i < 7.7 Gyr, to the 
(4) cooling time at this radius. Hence, M c ] ass j ca i(r) = M aas (r < 
'"cooi)/fcooi(?"cooi) is a measure of the rate at which the mass 
should be dropping out of the X-ray band provided there is no 
source of heating. As the main focus for this work is on the 
central parts of clusters where cooling is most dominant, the 
X-ray analysis was performed only for the core regions. 



3. Results 

3.1. CC and NCC cluster fractions: With and without a 
central radio source 

A well-known problem related to cooling-flows has been that 
of choosing an apt diagnostic for determining a cool-core clus- 
ter. To sep arate out the cool-co re clusters from the non-cool- 
core ones, iHudson et al.l d2008l) searched for a bimodality in 



several of the X-ray observed and derived quantities, such as, 
the central cooling time, f coo i, the cooling radius defined as 
the radius out to which f coo ] < 7.7 Gyr, the central entropy, 
the central density, the central luminosity, the mass deposi- 
tion rate, the central temperature drop and the slope in surface 
brightness profile, wherein "central" refers to 0.004/?5oo (0.4% 
of Rsoo). Based on the K-Mean Method (KMM) algorithm 
(lAshman et al.lll994l) as a test for bimodality (or tri-modality) 
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Fig. 5. Cooling time as an indicator of cool core. Left: The central temperature drop. Right: The slope in the surface brightness profile, 
cuspiness, at 4% of ^500- 



in the parameters, iHudson et alJ (I2008I) found f coo i as the best 
measure for cooling to divide the CC and NCC clusters. 

The f coo i-distribution histogram shown in Figure |4] shows a 
peculiarity in that the distribution is marked by a sharp increase 
at fcool < 1 Gyr; the fraction of clusters with f coo i < 1 Gyr be- 
ing 44%. This oddity at f coo i < 1 Gyr is also visible in two 
other properties, (1) the central temperature drop shown in the 
left panel of Figure [5] defined as the ratio of the temperature in 
the innermost region to the virial temperature and (2) the slope 
in the surface brightness profile, the 'cuspiness', shown in the 
right panel of Figure|5] Both the quantities show a break around 
a central cooling time of f coo i < 1 Gyr, even though the decrease 
in the central temperature drop is much more pronounced than 
the increase in the cus piness. This re s ult is a lso in concert with 
the recent findings of iRaffertv et al.l (12008b . who investigated 
the relation between star formation, cooling activity in the ICM 
and AGN heating, based on a sample of 47 cluster center domi- 
nant galaxies (CDGs). According to their study, only the CDGs 
with cooling times below 0.8 Gyr exhibit positive color gradi- 
ents, signifying an increase in star formation with decreasing 
distance from the galaxy center. The underlying reason for this 
behavior is not yet clear. It may be that the cool gas at the cen- 
ters of galaxy clusters is feeding the star formation, in which 
case the short cooling times should be tied with the time-scale 
over which the stars form and their light declines. It may well 
also be that AGN activity at the center of mass flows triggers 
star formation at the central regions. If the AGN feedback is re- 
sponsible for regulating the cooling flows in clusters, then the 
cooling time-scale may possibly reflect an intimate link with 
the time scale of the AGN outbursts (Section|4]i. 

Based on the above results, we divided our sample into 
three categories, (1) strong cool-core (SCC) clusters with 
fcooi < 1 Gyr, (2) weak cool-core (WCC) clusters with 1 Gyr< 
fcooi < 7.7 Gyr (the upper limit of 7.7 Gyr is the usually as- 
sumed value for the cooling time corresponding to z — 1, signi- 
fying the lookback time since the last major heating event, see 
McNamara & Nulsen 2007) and (3) non-cool-core (NCC) clus- 
ters with fgooi > 7.7 Gyr. These cuts result in 44% SCC clusters, 
28% WCC clusters and 28% NCC clusters. The need to divide 



the distribution into three subgroups is bolstered by the fact that 
the KMM test showed adding a third sub-group improved the 
likelihood ratio, giving rise to a tri-modal distribution. 

On cross-correlating the clusters with the presence of a 
CCRS, we find that all SCC clusters show cluster-center radio 
sources (Figure [5). The resulting fraction of CCRSs amongst 
WCC clusters is 67% and that in NCC clusters is 45%. A non- 
negligible fraction of CCRSs in the WCC and NCC cluster 
population makes it uncertain whether there is a fundamental 
one-to-one correspondence between AGN heating and the lack 
of the expected cluster cooling. On the other hand, the probabil- 
ity of a BCG manifesting AGN activity clearly increases with 
decreasing cooling time. The next question that then arises is 
whether the radio luminosity of the central cluster radio source 
itself is correlated with f coo i. This is shown in Figure [7] for 
SCC clusters and WCC clusters. From hereon, we refer to the 
combined set of SCC and WCC clusters as the cool-core (CC) 
clusters. This plot does not present a straight-forward interpre- 
tation of the interdependence between the AGN synchrotron 
power and the cooling time-scale. As a whole, there seems 
to be an anti-correlation between the two quantities but this 
seems to break down for clusters with t coo \ < 1 Gyr. This ap- 
parent anti-correlation could be indicative of a need for more 
powerful AGN as heating agents in clusters with short t coo \. 
Yet the absence of any correlation between the f coo i and Lr at 
short cooling times (< Gyr) implies that the AGN luminosity is 
more sensitive to a physical quantity other than the gas cooling 
time, such as possibly the mass deposition rate, M c iassicai (see 
Section [3.2b . 

Shown as crosses in Figure [7] are four systems, NGC4646, 
NGC1550, NGC5044 and MKW4, which clearly depart from 
this trend and all of which are groups. A general property that 
the groups in our sample seem to share is that apart from hav- 
ing low temperatures (both virial and central), they also all have 
high central densities and subsequently short f coo i (see Eq. 3). 
On the other hand, the groups tend to have very steep den- 
sity gradients resulting in small classical mass deposition rates, 
^classical- In other words, M c i ass j ca i is more sensitive to the mass 
encompassed within the integration radius (see Section 12.2. Il l 
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Fig. 6. The fraction of strong cool-core (SCC) clusters, weak 
cool-core (WCC) clusters and non-cool-core (NCC) clusters in the 
HIFLUGCS sample. Also shown are the fraction of clusters containing 
central radio sources for each category (shaded). 

than the cooling time at that radius [also see Figure 6(G) of 
Hudson et al. 2008]. The behavior of Lr versus M c iassicat is 
investigated in Section 13.21 The fourth outlier, MKW4, is an 
interesting cluster under intensive study at radio wavelengths 
(see Section 13.3. U . Assuming the anti-correlation interpreta- 
tion is correct, the best fit powerlaw excluding the four outliers 
derived using the bisector l inear regression routine, BCES from 
Akritas & Bershadvld 19961) is 



10 42 /j- 2 ergss- 1 



(0.041 + 0.016) x 



^cool 

Gyr 



-3.16±0.38 



(5) 



This routine, like FITEXY, includes uncertainties in both the 
quantities but also additionally performs the minimization in 
both the dimensions. The Spearman rank correlation coefficient 
is -0.63 and the probability for the null-hypothesis is 8 x 10 . 

For comparison with other works, we also determined the 
fraction of CCRSs in CC clust ers, the fraction being 8 7%. This 
is consistent with the result of Dunn & Fabian ( 20061) . who an- 
alyzed a low-redshift sample of clusters (B55) selected from 
pre-ROSAT data. Even though they find a slightly higher frac- 
tion (95%) of CC clusters with CCRSs, they used a lower cut 
in f coo i to determine CC clusters and, additionally, selected only 
those clusters which showed a central temperature drop > 2. 
Using these criteria reduces the fraction of CC clusters in our 
sample to 25% but incre ases the fracti on of CCRSs in CC clus- 
ters to 100%. Similarly. iBurns [1 19901) finds a somewhat lower 
fraction of 70% but the classification into CC and NCC clusters 
therein is based on the Hubble time. Using the Hubble time as 
the cut in f coo i increases the fraction of CC clusters in our sam- 
ple to 89% and reduces the fraction of CCRSs i n CC clusters 
to 78%. We also bear in mind that the result by IBurns (1990) 
is based on an incomplete sample a nd old X-ray Einstein data. 
Furthermore, the radio data used by BurnsI (Il990l) are based on 
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Fig. 7. The central cooling-time versus the integrated radio luminos- 
ity for the CCRSs in CC clusters (SCC + WCC). The black solid line 
represents the anti-correlation trend which breaks down for clusters 
with f coo i < 1 Gyr. The labeled clusters are outliers with peculiar prop- 
erties (see text for more). 



monochromatic 5 GHz VLA observations sensitive to largest 
structures of only about an arcminute, which in some cases 
might lead to over-resolved structures and, hence, an under- 
estimation of the radio luminosity. 

3.2. Cooling and AGN activity 

We looked for correlations between the radio luminosity of the 
CCRSs and the X-ray-derived quantities to allow us to identify 
the underlying mechanisms that link the AGN activity and the 
cooling properties in clusters. 

Shown in the left panel of Figure[8]is the bolometric X-ray 
cluster luminosity in the energy range 0. 01-40 keV as inferred 
from ROSAT and ASCA measurements dReiprich & Bohringei 
20021) . Lx, versus the integrated radio luminosity (see 
Section I2T2I for the 48 clusters with CCRSs. For the SCC 
clusters, shown as filled (blue) circles, there is a clear positive 
trend visible, although with a considerable spread. Since the 
X-ray luminosity is related through scaling relations to other 
global parameters of a cluster, such as the T vn and cluster mass, 
similar correlations may be obtained between the radio power 
of a CCRS and these quantities. This is the first time that the 
radio power of a centrally located AGN, the prime candidate 
for counteracting the cooling of the X-ray radiating ICM gas, 
has been shown to be correlated with the large-scale cluster 
properties. This result implies that there is a link between two 
regions, vastly differing in scales; the region over which AGN 
accretion takes place, which is no more than a few hundredth of 
a parsec, and the ICM, which extends out to 1 - 2 megaparsecs. 
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Fig. 8. Radio and X-ray correlation plots. Left: Total radio luminosity vs. ROSAT bolometric X-ray luminosity for SCC (filled blue circles), 
WCC (open green circles) and NCC (open red squares) clusters. Right: Total radio luminosity vs. the classical mass deposition rate for SCC 
clusters (filled blue circles) and WCC clusters (open green circles). The black triangles are four WCC clusters with M d;lssic;l i > 1 and no central 
radio source. Hence, these are only upper limits. The filled red square on top left is A2634, also a WCC cluster. The solid black line is the best 
fit through CC clusters (SCC and WCC clusters combined) excluding A2634. 



Using the two-dimensional bisector linear regression routine 
BCES, the trend between the Lr and Lx for SCC clusters may 
be quantified as below: 



lO 42 /!^ 2 ergss- 



a x 



10 44 h~\ ergs s" 1 



(6) 



where a = 0.03+0.01 andb = 1.38 ±0.16. The Spearman rank 
correlation coefficient of the fit is 0.64 and the probability of 
their being no correlation is 2.7 x 10~ 4 . The Pearson correlation 
coefficient for the fit is larger and is equal to ~ 0.76. 

Note that luminosity-luminosity plots should be considered 
with caution d ue to the common redshift-dependence in both 
the quantities dKembhavi et al.lfl986t lAkritas & Siebertlll996t 
Merloni et alj 120061) . Even though there are no censored data 
points (upper limits) for the category of SCC clusters, in that 
every SCC has radio source at the center, spurious correlation 
may still be introduced due to the common dependence on the 
distance. In order to check for such an occurrence, we simu- 
lated randomized radio and X-ray luminosities confined to the 
observed ranges follow ing the distributions, n(Lx) dLx °c L y ° 7 
dBohringer et all l2002h and n(L R )dL R oc L R 78 dNagar et al 



20051) . where n is the source number density. We assigned ran 



domly distributed redshifts to the randomized luminosity data 
sets, according to the law n ~ D^. These luminosities were re- 
observed after applying the X-ray flux limit, fx (0.1-2.4) keV> 
2 x 10" 11 erg s _1 ctrT 2 , the same as that for the HIFLUGCS 
sample, and the radio flux limit, 1.5 mJy the average point- 
source detection limit for NVSS. The resulting Pearson cor- 
relation coefficients, p P , were compared to the observed one. 



Based on these simulations, we compute the probability of hav- 
ing pp > 0.76 and the correlation slope equal to or greater than 
that observed for a completely randomized set of X-ray and 
radio luminosities as less than 1% (a spuriously induced cor- 
relation should produce a slope of around unity). This fraction 
increases to ~ 2.5% if instead the Spearman correlation coef- 
ficient is used and, if in addition, the observed probability of 
null hypothesis is used as a further constraint, i.e., the proba- 
bility of null hypothesis for the simulated data sets should be 
lower than 2.8 x 10~ 4 , then this fraction decreases to ~ 1.5%. 
Thus, we conclude that the probability of the observed correla- 
tion between the cluster X-ray luminosity and the radio lumi- 
nosity of a CCRS to be spurious is very unlikely. However, that 
such an induced correlation is possible at a level of ~ 3% in 
the worst case scenario, is worth keeping in mind for past and 
future studies on similar topics. 

Shown in the right panel of Figure[8]is the radio luminosity 
versus M c iassicai- This plot shows an even stronger trend than 
that seen with Lx- This further strengthens the likelihood of a 
coupling between gas cooling and the magnitude of the AGN 
activity. The NCC clusters do not appear on this plot since these 
by definition have no cooling radius, that is the central cooling 
time for these clusters is greater than 7.7 Gyr, implying zero 
mass deposition rates. 

There are two interesting subsets of clusters pertaining 
to the right panel of Figure [8] which deserve attention. The 
first subset comprises clusters which lack a CCRS but have 
^classical > 1 M yr _1 , and the second subset, not shown 
in Figure [8] comprises clusters which have a CCRS but for 
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which M c i ass i ca i = 0. The former subset (denoted by black tri- 
angles in Figure IS) consists of: A1650, A2589, A2657 and 
A 1060, with mass deposition rates (93.7 + 28.2) M Q yr -1 , 
(19.9 + 12.3) M yr\ (15.4 + 10.9) M yr~ ! and 
(6.6 + 3.O)M yr~' respectively. In order to understand the be- 
havior of the above four clusters, which are all WCC clusters, 
and to find out whether there exists a quantity that separates 
them from the other WCC clusters with a central radio source, 
we examined the entropy profiles of these clusters. Entropy is 
a powerful tool which provides information about two cluster 
parameters simultaneously - the temperature and the density; 
K(r) = kT(r)n{ry 2 ^, where r is the radius from the clus- 
ter center. Shown in Figure [9] are the entropy files of all but 
two WCC clusters. The two exceptions are A3266 and A3667, 
which have no CCRS but also have classical mass deposition 
rates consistent with zero. In other words, these are cases at the 
border line between WCC and NCC clusters and have, there- 
fore, been excluded from Figure [9] As can be seen, the entropy 
profiles of these four clusters (shown as color curves with sym- 
bols other than '+' symbols) are in no sense different from the 
rest. If anything, the profiles of A1650, A2589 and A1060 seem 
to continue to fall with decreasing clustercentric distance. This 
implies a steady increase in the gas density with decreasing 
radius in these clusters and, hence, relatively strong cooling. 
What is the source of heating in these clusters? 



Donahue et al.l (120051) investigated one of the above radio- 
quiet CC clusters, A 1650, using Chandra observations and 
proffered one of the following two explanations for the ab- 
sence of a radio AGN at the cluster center; (1) either the cluster 
has not reached the point where heating is necessary, or (2) the 
cluster experienced a major heating event about 1 Gyr ago such 
that it has not required feedback since then. Their conclusions 
are based on a lack of central temperature gradient in A 1650 
and a markedly raised central entropy as compared to other 
CC clusters with radio emission. Although the Chandra ob- 
servations used by us (including ~ 200 ks that became publicly 
available in 2008) also imply an insignificant central tempera- 
ture drop (7o/r v j r ~ 0.8), the estimated central entropy is not 
any higher than the average central entropy of the rest of the 
WCC clusters (Figure IS . 

Although, all of the above four clusters pose a serious 
threat to the AGN-regulated feedback fabric in cool-core clus- 
ters, A 1650 is most intriguing due to a high value of discrep- 
ancy between the expected and measured mass deposition rates 
[Medical ~ (93.7 + 28.2) M yr 1 and M spec < 0.7 M yr" 1 , 
where M sp Jlis the spectral mass deposition rate]. Interestingly, 
that there has been a mention of a weak detection of a ra- 



dio source at the center of A 1650 by Dunn & Fabianl (120061) . 
which in turn is base d on the VLA detection at 327 MHz by 



Mar kovic et al.l (120041) . who give the total flux-density of the ra- 



dio source at this frequency as 59 mJy. But we have been unable 
to re-confirm th i s claim using the same observations as used by 



Markovic et al.l (2004) down to 3 mJy , three times the back 



ground noise. As also pointed out by iDonahue et al.l (120051) . 
there is neither an indication of a past AGN outburst, either 




Fig. 9. The entropy profiles of the WCC clusters. All with a central 
radio source are shown in black (plus) symbols and four without are 
shown in colored (non-plus) symbols. The errorbars on those with a 
central radio source are omitted for clarity. The strong jumps apparent 
in some of the entropy profiles are insignificant relative to the error- 
bars. 



in the form of low-brightness diffuse lobe emission or cavities 
in the X-ray emission, nor are there any signatures of a recent 
merger. This cluster, along with its three companions, deserves 
further study in order to analyze other possible sources of heat- 
ing such as conduction, intracluster supernovae or preheating. 

The second subset corresponding to clusters with zero clas- 
sical mass deposition rates but which contain a CCRS consists 
of: A3391, A3395s, A3376, A0400, A1656, A3158, A2147, 
MKW8 (in decreasing order of L R ). These are NCC by defi- 
nition and show signs of cluster mergers at diffe rent stages, as 



do mo st of the other NCC clusters in our sample dHudson et al 



2008), based either on their X-ray properties (presence of sub- 
clumps or non-negligible separation between the X-ray emis- 
sion peak and emission weighted centroid) or radio proper- 
ties (presence of radio halos or relics or both). But contrary 
to the clusters in this subset, the remaining NCC clusters are 
devoid of a central radio source, supporting the idea that the 
presence of AGN is tightly correlated to gas cooling in clus- 
ters. The obvious question that then surfaces is, how does 
this subset fit into the AGN-heating and gas-cooling machin- 
ery? There are several solutions to this apparent contradiction. 
Firstly, a well-known fact - AGN are found at all locations in 
clusters and not only at the cluster centers (although with an 
inc reasing probability with decreasing cluster centric distance, 



2 A detailed descr iption of how M spC c is calculated can be found in 
iHudson et all J2008h 



see iMorrison & Ow en 2003; IBest et al.ll2007l) . There also ex- 
ist field- AGN with no apparent reservoir of bulk cool gas, such 
as is available in clusters. Hence, there evidently are mecha- 
nisms other than those related to cluster cooling that can trigger 
radio nuclear activity in galaxies. Secondly, there is evidence 
that mergers may play a role in activating the central engines 
of AGN by transferring gas to the cluster galaxies and provid- 
ing material for bo th, AGN accretion and also star formation 
(lOwen et al. I ll999h . But alongside there also are contradictory 
findings according to which mergers may as well strip away 
gas from galaxies and result in inhibition of both the processes 
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(e.g. see lGiacintucci et al 1 l2004l) . Thirdly, a configuration con- 
taining a NCC cluster with zero mass deposition rate and a cen- 
tral radio source may be obtained through a merger between a 
NCC cluster and a CC cluster, latter harboring a central AGN, 
in such a way that it results in disruption of the cool core and 
leaves behi nd only trac e s of th e past cooling activity. The sim- 
ulations of iBurns et alj d2008l) indeed show that NCC clusters 
are formed as a result of major mergers right at the beginning of 
their evolution, whereby they grow in time at the expense of CC 
clusters. As an example, A2634 is a WCC cluster based on the 
central cooling time but in most other respects is closer to being 



a NCC cluster than a WCC cluster. According to lHudson et al 



d2008l) . the X-ray morphology of this cluster is consistent with 
that of a merging cluster. Yet it has a cool core (short f coo i), and 
so A2634 might be a strong candidate for such a case where 
the cool core is being destroyed by a major merger. 

Excluding the cluster on the upper-left corner (A2634) of 
Figure [8] and the aforementioned subsets of outliers, the pow- 
erlaw fit for SCC and WCC clusters using the BCES routine 
gives, 



10 42 /^ 2 ergss- 



/ ^classical \ 

= a x - 

I Mo yr-V 



(7) 



where a = (4.7 + 5.0) x 1(T 5 and b = 1.69 ± 0.25. It can be 
argued that this exercise may be more appropriate if we use 
the spectrally determined mass deposition rate, M spec , as op- 
posed to the classical one, M c lassical> as the former gives the 
actual observed rate at which the gas is cooling out and accret- 
ing onto the supermassive black hole. M spec is, in fact, the fuel 
for the central AGN and should be correlated with the AGN 
output. On the other hand, the question that we are trying to 
address here is whether the AGN output can account for the 
difference between the classical (predicted) and spectral (ob- 
served) mass deposition rate and, thereby, provide a solution 
to the cooling flow puzzle. In addition, excluding clusters for 
which M spec is consistent with zero, the classical mass deposi- 
tion rates exceed the spectral mass deposition rates on average 
by a factor of 15. Hence, (M c i assica i - M spec ) ~ M c i assica i and the 
correlation between the radio luminosity and the two types of 
mass deposition rates is expected to be similar. We also point 
out that the spectral mass deposition rates estimated using the 
ACIS instrument on Chandra are, due to its low spectral re- 
solving power, only moderately accurate. Therefore, any cross- 
correlation with the spectral mass deposition rates will have a 
large uncertainty making robust interpretations difficult. 

Finally, we interpret the strong correlation seen between 
^classical and Lr as supporting evidence for a feedback system 
in which the AGN activity is more enhanced in clusters with 
higher mass deposition rates, and the AGN in turn quenches 
the cooling of gas by heating the am bient medium. We note 
that a similar result was obtained by iPeres et alj d 1 998b be- 
tween 5 GHz spectral radio luminosity and the mass deposi- 
tion rate determined with ROSAT, although for a much smaller 
subset constitutin g only 15 galaxy clusters of the B55 sample 
dEdge et alJll990h . 



3.3. Brightest cluster galaxies (BCGs) 

The brightest cluster galaxies are unique in terms of their high- 
luminosity and proximity to the centers of their host clusters. 
The BCGs are extremely interesting objects and have long been 
subjects of a wide range of studies. At one extreme, their for- 
mation and evolution is closely tied with the Mpc scale clus- 
ter environment in which they reside. At the other extreme in 
the hierarchy of structure formation, BCGs are just one level 
above the ~ 10 -4 pc scale supermassive black holes (SMBH). 
The BCG bulge properties, such as the optical bulge luminos- 
ity and stellar velocity dispersion, obey certain scaling relations 
that permit indirect estimation of the mass of SMBH s (e.g. 



Kormendv & Richstonel 19951 iFerrarese & MerritthoOOl) . 

In this section, we correlate the BCG magnitudes or, equiv- 
alently, the mass of the SMBHs with the AGN radio luminos- 
ity and the large-scale X-ray properties of galaxy clusters. The 
mass of the SMBH, Mbh, is derived using the scaling relation 
between the near-infrared ( NIR) bulge magn it ude a nd the in- 
ferred Mbh, as deduced by iMarconi & Hunt! (120031) . Further, 
we test whether there are any distinctions in the BCGs proper- 
ties amongst the three different types of clusters (SCC, WCC 
and NCC clusters). 

The BCG apparent magnitudes were taken from the 
Two Micron All-Sky Survey (2MASS ISkrutskie etai1l2006l) 



Extended Source Catalog (XSC0). We used the 2MASS total 
magnitudes in /if -band (2.16 /urn), km ext , estimated from ex- 
trapolation of the surface brightness profiles (SBP). The ex- 
tended source detection limit for /if -band at 10-cr is 13.5 mag 
and the uncertainties range from 0.02-0.23 mag with a mean 
of 0.06 mag. In short, the 2MASS SBPs have been derived 
from fitting a modified Sersic function to the elliptical radial 
light distribution of the BCGs. The total magnitudes are es- 
timated from summing two terms. The first term corresponds 
to the isophotal magnitude estimated from fitting an ellipse 
to the standard isophote of mean surface brightness of K20 = 
20 mag arcsec~ 2 . The second term is derived by integrating the 
best-fitting Sersic law starting from the standard isophote, r2o, 
out to a delimiting isophote, which is typically about four scale 
lengths. This 2MASS strategy insures that the total flux of a 
galaxy is recovered. 

The BCGs in 2MASS were located by searching around 
the brightest cluster galaxy in each cluster using an initial com- 
pilation, kindly provided by Heinz Andernach, based on the 
available data at NASA/IPAC Extragalactic Database (NEE0) 
and Hyperleda5 All the BCGs were found within 5 of the 
given search position except in the case of three clusters, 
A2204 (SCC), A2065 (WCC) and A2163 (NCC), where, based 
on visual inspection, 2MASS did not manage to locate the 
right galaxies (i.e. even though within 5 of the given search 
position, the putative BCGs were not found in the 2MASS- 
XSC catalog). For these three clusters, we retrieved the 2MASS 
/if -band Atlas images which have a plate scale of 1 /pix, 
and fitted the BCG surface brightness distribution s using the 
two-dimensional galaxy-fitting program, GALFIT dPeng et al 



3 http://irsa.ipac.caltech.edu/Missions/2mass.html 

4 http://nedwww.ipac.caltech.edu 



http://leda-univ-lyon 1 .fr 



12 



Rupal Mittal et al.: AGN heating in the HIFLUGCS sample of Galaxy Clusters 



20021) . In order to be consistent with the 2MASS fitting rou- 
tines, we restricted the functional form to Sersic models. While 
the SBPs of the BCGS of A2065 and A2163 were well- 
reproduced with a single-component Sersic models, the BCG 
in A2204 required a double-component Sersic model. For 
comparisons of the /if -band magnitudes between the 2MASS- 
XSC and GALFIT estimates, we carried out tests by apply- 
ing GALFIT to a few of the BCGs present in the XSC catalog 
and, hence, with known 2MASS magnitudes. We found that 
whereas GALFIT systematically underestimates the magnitude 
for bright BCGs (ktn ext < 10 mag) by about 10%, the GALFIT 
magnitudes are consistent with the 2MASS-XSC magnitudes 
for faint BCGs {km ext > 10 mag) to within 3%. Since the 
GALFIT magnitudes of the BCGs in A2204, A2065 and A2163 
are all fainter than 12 mag (but brighter than the detection limit 
of 13.5 mag), we deem them to be trustworthy to within 5%. 

The apparent magnitudes from 2M ASS were corrected fo r 
Galactic extinction using values from ISchlegel et al.l d 1 998b . 
the typical correction values being small - on the order of 
~ 10 2 , and were then converted into absolute magnitudes us- 
ing the redshifts compiled from NED. We did not apply any k— 
correction since these galaxies are all nearby. 

3.3.1 . Supermassive black hole mass and radio 
luminosity of the BCG 

There has been a lot of debate over the use of BCG scaling 
relations for determinations of Mbh- This is due to the fact 
that many BCGs are often accompanied by low surface bright- 
ness env elopes extending out t o, as far as, several hundred kilo- 
parsecs (iGonzalez et al]|2005l) . These are the well-known "cD 



galaxies". The extended envelope s, also known as intraclus- 
ter light (ICL) dLin & Mohrj|2004 . are thought to either rep- 
resent debris accumulated over the merger history of the BCG 
or from tidal stripping of other cluster galaxies. The extended 
emission may as well originate from stars forming out of the 
condensed gas in cooling clusters. An important investigation 
we will follow-up in a subsequent paper is to study how the 
luminosity of these envelopes correlates with cooling parame- 
ters, such as the f coo i, M c i ass ; ca i, M spec etc. This requires careful 
decomposition of the BCG light profile into an inner compo- 
nent, associated just with th e galaxy, and an out er flatter com- 
ponent representing the ICL. Seigar et aU ( 2007 ) fitted analyti- 
cal models with two Sersic components to separately measure 
the profiles of the central and extended parts of 5 cD galaxies 
and showed that the contribution of the envelopes to the total 
light is around 60% to 80%. 



De spite the above factors, recent studies ( IBatcheldor et al 



2007); iMarconi & Huntl I2003I) have shown that the SMBH 



masses derived from the scaling relation using the NIR bulge 
magnitudes of the BCGs are consistent with those derived 
from stellar velocity dispersions (cr v ), both the relations yield- 
ing similar amount of scatter. This is not quite true for V- or 
Z?-band magnitudes though, where the scatte r is much higher 



compa red to the preferred <x v -Mbh relation. IBatcheldor et al 



(I2007I) attribute the scatter to inclusion of luminosity from the 
outer-envelopes of the cD galaxies, which may be more pro- 
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Fig. 10. The mass of the SMBH versus the total radio luminosity of 
the BCG in SCC (blue filled circles), WCC (green open circles) and 
NCC (red open squares) clusters. Visible is a weak trend given by solid 
blue line for the SCC clusters such that the Lr of the BCG increases 
with the SMBH mass. Labeled are the five SCC outliers. 



nounced in the V- or B- band, especially if the outer-envelope 
light represents on-going star-formation, and which might not 
ha ve anything to do w ith the central galaxy dynamics (however, 
see lLauer et al.ll2007l for a difference of opinion). Based on the 
studies on the NIR bulge magnitude-MBH relatio n, we use the 
following scaling relation ( Marconi & Huntl2003 ) to derive the 
SMBH masses : 



log 



Ml 



Mr 



M n 



log 



10 ~i 



10.9 



(8) 



where a = 8.21 + 0.07 and b = 1.13 + 0.12. In order to convert 
the absolute magnitudes into luminosities in units of /f-band 
solar luminosity, we used the absolute ^T-ba nd solar magnitude 
equal to 3.32 mag dColina & Bohiinlll997l) . In Figure [TOl we 
present the SMBH masses versus the Lr for the 48 CCRSs 
(BCGs) in our sample. Whereas on the whole there appears 
to be a poor correlation of increasing AGN radio output with 
increasing Mbh, categorization of clusters based on f coo i results 
in a trend to appear between the two quantities but only for the 
SCC clusters. The best-fit powerlaw using the two-dimensional 
fitting algorithm BCES for the SCC clusters is: 



I0 42 h- 2 ergss- 



(0.008 + 0.004) x 



M] 



BH 



10 9 Mr, 



.10±0.42 



(9) 



The Spearman rank correlation coefficient of the fit is 0.46 
and the probability of their being no correlation is 1.3 x 10~ 2 . 
The Pearson correlation analysis yields a larger correlation co- 
efficient (0.59) and a much smaller probability of the null- 
hypothesis (9 x 10~ 4 ). 
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Fig. 11. The 2MASS K-band total magnitudes of the 64 BCGs in the HIFLUGCS sample versus M 500 in the left panel and ROSAT bolometric 
cluster X-ray luminosity in the right panel. The different symbols and colors denote SCC (filled blue circles) clusters, WCC (non-filled green 
circles) clusters and NCC (non-filled red squares) clusters. The solid (black), the dashed (blue) and the dotted (pink) lines denote the best-fit 
powerlaws for all clusters, only SCC clusters, and only non-SCC (WCCs+NCCs) respectively. 



Several studies have explored the correlation between the 
SMBH mass and radio luminosity of the embedded AGN, lead- 
ing to contradictory results. Whereas it is clear that the radio 
loudness, the ratio of the radio to the optical lumi nosity of an 
AGN, is a strong function of the black hole mass ( ILao: 



try ot a n 

ji feoori . 



whether or not there exists a tight relation between L r and Mbh 
is still a matter of debate. Franceschini et al.l d 1998b derived a 
tight relation between the radio luminosity at 5 GHz, L5 qhz, 
and Mbh using a sample of 13 nearby early-type weakl y ac- 



tive galaxies, such that L5 ghz k tfff. lLacy et al.l d2()()l ) also 
find a similar correlation for a sample of steep-spectrum ra- 
dio sources , albeit leading to a flatter slope of 1.4. Although 
Laori d2000l) confirmed this trend of increasing radio luminos- 
ity with increasing black hole mass using a much larger sample 
of 29 nearby galaxies and 89 PG (Palomar-Green) quasars, the 
trend they obtain is weak an d presents a lot of scatter. As an 
example of dissenting views, iLiu et al. I d2006l) find no relation 
of L5 ghz against the black hole mass. On the other hand, they 
find a strong correlation between the jet power and the black 
hole mass, where they used the flux density of the radio lobes 
at 151 MHz as a proxy to determine the jet power. Based on 
the above results it is not yet clear how the radio luminosity of 
an AGN scales with the black hole mass. It may be that differ- 
ent black hole efficiencies, conversion rates from the total AGN 
power to the radiative power of the jet and large-scale environ- 
mental effects from source to source, cause the radio luminosity 
to display a wide range of values for a given black hole mass. 

In pursuit of determining whether or not there is an under- 
lying physical mechanism that ties the radio power to the mass 
of the SMBH, Figure [10] presents an interesting outlook. We 
investigate for the first time the dependence of Lr on Mbh tak- 
ing into consideration the different environments (SCC, WCC 
and NCC clusters) in which the central ly located radio sou rces 
reside. According to the simulations of lBurns et al.l (120081) . CC 



clusters grow an enhanced cool-core gradually and steadily via 
small mergers, unlike NCC clusters which experience major 
mergers early in their history. Hence, repetitive small mergers 
in addition to providing material for the growth of the cool core 
in CC clusters, may also result in a quasi-steady mass accre- 
tion rate onto the SMBH in CC clusters. The fact that only the 
CCRSs in SCC clusters show some level of correlation between 
Mbh and Lr implies that the AGN output in these systems is 
proportional to the average mass-accretion rate onto the black 
hole, thereby, balancing the radiative losses. 

At the time of writing this article, MKW4, an outlier below 
the best-fit line in Figure[lOl had no low-frequency (< 1 .4 GHz) 
radio data available. This inspired us to acquire 327 MHz VIA 
data for this cluster, which is work-in-progress. MKW4 is an 
interesting cluster in that high-frequency (1.4 GHz and 5 GHz) 
VIA archival data showed a cluster of closely spaced point- 
sources > 1 arcmin to the NE of the BCG radio emission. 
Furthermore, radio emission from the BCG was only detected 
at 5 GHz and not at 1.4 GHz, the latter having had only 2- 
minute on-source integration time. The preliminary 327 MHz 
image of MKW4 revealed the same feature NE of the BCG but 
due to insufficient resolution, it still remained unclear whether 
this emission was associated with the nucleus of the CCRS or 
corresponded to a high-redshift cluster system; the latter con- 
clusion being based on the proximity of these point sources 
to each other. Assuming that the emission in the NE reflected 
the region of intense interaction between the lobe and ICM and, 
hence, was a part of the central radio source, would have caused 
the Lr to increase 30 fold and for it to no longer be an outlier. 
To confirm this hypothesis, we acquired dynamic VIA time at 
1 .4 GHz in CnD configuration with 3 hrs integration time. Even 
though the new 1 .4 GHz data show no signature of a connec- 
tion between the CCRS and the bunch point-source emission in 
the NE, it clearly shows emission from the CCRS, which went 
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undetected with the old archival data. This itself increased the 
total radio output by more than an order of magnitude bringing 
this cluster closer towards the observed M BH -L R trend for the 
SCC clusters (this change is incorporated in Figure ITOb. The 
follow-up study of the CCRS in MKW4 underlines the impor- 
tance of obtaining reliable spectra of radio sources. 

In addition to MKW4, there are four other outliers in 
Figure [10] which lie above the best-fit relation: A0780, A2597, 
MKW3S and A3581. We believe these systems might be ones 
which have experienced powerful radio outbursts in recent past 
resulting in the present AGN heating rate to be greater than 
the average rate at which the mass has been accreting onto the 
SMBH since the formation of the cool core. Similarly, MKW4 
may reflect a CCRS that has been caught just at the beginning 
of another heating cycle and has yet to reach its peak radio ac- 
tivity. These outliers, both below and above the best-fit line may 
be reflective of the episodic nature of CCRSs in some systems. 
Another plausible reason for these outliers may be the weak- 
ness in the underlying assumption of the integrated luminosity 
of the BCG as a robust indicator of the bulge luminosity in all 
the cases. 



3.3.2. Large-scale cluster properties and BCGs 

Inherent searings between BCGs and clusters h ave been im 



clusters. A possible reason for such a separation may be re- 
lated to the continuous growth of BCGs, as already mentioned 
in Section D. 3. II due to iterative small-scale mergers through 
which the CC clusters grow. The BCGs in NCCs, on the other 
hand, form constituents of non-relaxed cluster environments 
with a history of one or more major mergers, and the subse- 
quent heating at the central regions may hinder further BCG 
growth at the same rate as that of BCGs in CC clusters. This 
argument is supported by the fact that on comparing the ra- 
dial profiles of nu merically simulated CC and NCC clusters, 
Burns et al.l d2008l) observe an excess of baryons in CC clusters 
relative to NCC clusters. 

An interesting note is that this segregation appears only 
between Lbcg versus M500 and not between Lbcg versus Lx. 
However, a similar intrinsic separation for CC and NCC clus- 
ters is seen between T wn and Lx- This is attributed to the fact 
that at a given temperature, SCC clusters have a higher lumi- 
nosity as compared to non-SCC clusters due to an increased 
gas density at the center. This causes the SCC clusters to form 
an envelope towards the higher luminosity end. Since M500 has 
been derived from T v [ r , it may be that the resulting magnitude 
of separation between M 5 oo and L x cancels with that between 
Lbcg and M 500 . 

The estimated best-fit powerlaw for the Lbcg-A^500 relation 
based on the BCES algorithm is: 



plied in numerous obse rvational studies (e.g. iLin & M ohr 



2004tlBroughet al 



(Zhen g et alj|2007 



ugh 

HE 



2008 ) as well as cosmologica l simulations Lbcg 



Coorav & Milosavlievicll2005l) . In this con- 10 1 1 hz} L 



= a x 



Ms 



10 14 M yr- 



(10) 



text, formation and evolution of BCGs and its dependence 
on the host cluster is an important tool to understand these 
scalings observed between the BCG luminosity and the host 
halo mass and X-ray luminosity. There are several proposed 
BCG evolution scenarios to support these observations, such 
as (1) dynamical-friction governed galactic cannibalism, (2) 
rapid mergers between galaxies during the epoch of cluster for- 
mation, (3) co-evolution of BCGs with cluster growth due to 
mergers embracing the paradigm of hierarchical structure for- 
mation and (4) cooling-flows. While the first two scenarios may 
cont ribute significantly d uring the early epo c h of B CG forma- 



tion, ILin & Mohri (120041) and iBrough et al.l (120081) argue that 
the BCGs co-evolve with the host clusters via mergers with 
the BCGs of the falling subclusters, which lead to subsequent 
growth of the BCG luminosity with increasing cluster mass. 

We present in this section the BCG-host cluster correlations 
for the HIFLUGCS sample, with the aim to extend previous 
studies by analysing a low-z flux-limited sample unique in its 
completeness and homogeneity in the ways of obtaining the X- 
ray and NIR quantities. The following study is unique in an 
additional aspect in that we bear in mind the possibility of dif- 
ferent growth histories for the BCGs corresponding to CC and 
NCC clusters. 

Shown in Figure QT| is M500 versus Lbcg in the left panel 
and Lx versus Lbcg in the right panel. Shown in different sym- 
bols and colors are the three different types of clusters. Whereas 
both the panels clearly indicate that the BCG grows with the 
cluster size, the left panel additionally shows a segregation be- 
tween the SCC clusters (blue filled circles) and the non-SCC 



where a = 3.525 + 0.277 and b = 0.624 + 0.054 for all clusters, 
a = 4.305 + 0.290 and b = 0.616 ± 0.005 for the SCC clusters 
only and a = 2.552 ± 0.362 and b = 0.752 ± 0.095 for the non- 
SCC clusters (CC clusters) only. The Spearman rank correla- 
tion coefficients for the fits are 0.67, 0.87 and 0.62 for all, SCC 
and non-SCC clusters, respectively. Looking at the fit results 
for all clusters and SCC clusters only, it is clearly seen that the 
SCC clusters have a higher normalization by about 20%(since 
the two fits have the same slope within the l-cr errorbars, the 
normalizations may be directly compared). The slope of the 
obtained Lbcg-A^soo relation is steeper than the values derived 
in other works, which tend to center arou nd 0.3, even though 
presenting a wide range from 0.1 to 0.5 dLin & Mohr! 120041 
Brough et al1l2008l:|Popesso et al.ll2007l:lHansen et al.ll2007l) . 



Going a step further, we may combine Equations (0 and 
(ITOb to derive a relation between the SMBH mass and the clus- 
ter mass: 



M, 



BH 



10 9 M, 



= a x 



M500 



10 14 M yr- 



(11) 



where a = 0.98 + 0.08 and b = 0.61 ± 0.06 for all clusters, a = 
1.15+0.09 and b = 0.63+0.06 for the SCC clusters only and a = 
0.75 ±0.14 and b = 0.70 ±0.12 for the non-SCC clusters. Even 
though indirectly derived using the BCG bulge luminosities, 
such a correlation could be indicative of a fundamental relation 
between the host cluster halo and central SMBH similar to the 
galaxy bulge mass-black hole mass relation. 
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Similarly, the best-fit powerlaw for the Lbcg-Lx relation is 
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= a x 



lO 44 ^" 2 ergss" 1 



(12) 



where a = 4.54 + 0.34 and b = 0.36 + 0.03 for all clusters, 
a = 5.15 ± 0.38 and b = 0.32 + 0.03 for the SCC clusters only 
and a = 3.49 ± 5.09 and b = 0.50 + 0.07 for the non-SCC 
clusters. The Spearman rank correlation coefficients for the fits 
are 0.68, 0.82 and 0.55 for all, SCC and the non-SCC clusters, 
respectively. The Lbcg-^x and Lbcg-^soo relations show that 
there is a statistically significant difference in slopes between 
SCC and non-SCC clusters. Similar correlations were obtained 



by iKatavama et al.l (12003b using a larger sample, which also 
contains the HIFLUGCS sample. However, their results show 
much higher scatter between the BCG luminosity-host halo 
mass, which is due likely to the use of optical B-band mag- 
nitudes which have larger error bars of mean value ~ 0.2 mag. 
Further, IKatavama et al.l 02003b correlated the optical magni- 
tude versus the total mass of the host cluster defined as the in- 
tegrated mass within a fixed metric radius of 5 Mpc, for which 
both, statistical as well as systematic uncertainties are expected 
to be larger. 

The above results, which highlight a strong dependence 
of the BCG NIR magnitudes on the scale of th e host cluster 



howev er, are at variance with those obtained by IBrough et al 
d2002l) . They claim that any division between BCGs in low- 
Lx and high-Lx clusters as seen in high-z clusters, disappears 
for cluster s with z < . 1. The input low-z sample (z < 0.1) 



studied bv lBrough et al.l d2002l) consisted of 150 Abell clusters 
with a flux limit of f x (0.1 - 2.4) keV= 3 x 10~ 12 erg s" 1 cnr 2 
in the ROSAT hard band (0.5-2.0 keV). After matching these 
with the 2MASS catalog resulted in a final sample comprising 
76 clusters with only those BCGs which have robust 2MASS 
magnitude measurements in the /iT-band. We argue that the con- 
tradiction in results might stem from differences in the aper- 
ture radii used, within which the magnitudes are calculated. 
The 2MASS database pro vides galaxy magnitu des based on 
a suite of apertures. While IBrough et al.l 00021) employ inte- 
grated magnitudes measured using circular apertures of a fixed 
metric radius of 12.5 hz\ kpc, we use, as explained above, 
the total aperture radii by extrapolating the SBPs, the mean of 
which for the BCGs in our study is ~ 50 }c\ kpc. From this 
we conclude that the total BCG magnitude is a better quantity 
to use, as opposed to the BCG core magnitude, for detecting 
correlations with global cluster properties. 



4. Discussion 

The riddle of cooling-flows in clusters has continued to baf- 
fle us. Recent high-resolution Chandra images revealed radio- 
loud AGN embedded in the centers of cool-core (CC) clusters 
surrounded by regions emptied of the X-ray emitting gas, sug- 
gesting a strong tie between the cluster central radio source 
and the cooling of ICM. Since the discovery of numerous 
AGN -blown bubbles in the atmospheres of CC clusters, various 
modes of energy transfer from the AGN to the ICM have been 
investigated. While the most successful mode has turned out to 



be dissipation of energy stored in the radio bubbles as they rise 
buoyantly through the ICM, previous studies have fallen short 
of a thorough investigation of the relation between the radio lu- 
minosity of the centrally located AGN and cooling properties 
of a CC cluster. 

It is now a widely accepted fact that the AGN activity is 
triggered by gas accretion onto the central black hole. Playing 
devil's advocate, it may then be argued that the AGN output 
is only to be expected to scale with the mass accretion rate. 
Under the assumption that the cool gas flows from the outer 
cluster regions to the very centers of the BCGs and serves as 
the fuel for the black hole, it is not surprising that the radio 
luminosity of an AGN should scale with the cluster mass, and 
also, even though to a lesser degree, the inverse of the cool- 
ing time of the gas. This explains the underlying trend seen in 
Figure [8] between the M c i ass ; ca i and the radio output of CCRSs 
in CC clusters. However, the picture thus developed so far does 
not contain any ingredients reflecting on a self-regulated cycle 
formed between gas cooling, star formation and AGN heating. 
In other words, it may well be that even though AGN activity 
enhances with cluster scale, the cooling of ICM is regulated by 
an altogether different process, such as cluster mergers. 

The first strong argument in favor of AGN-regulated heat- 
ing comes from the observation that the AGN fraction increases 
with decreasing central cooling time, f coo i being the best di- 
agnostic to distingui sh CC from NCC clusters. The study by 



Raff ertv et al.l d2008l) shows that the central star-formation rate 



also is a strong function of f coo i (see Section 13. 11 1. That only 
clusters with short f coo i < 1 Gyr have increasing on-going 
star formation with decreasing clustercentric distance, implies 
a chain of intricately linked processes which maintain heat- 
ing and cooling rates in cluster atmospheres. These results to- 
gether call for a feedback process in which AGN heating be- 
comes more of a requisite in clusters with shorter cooling times. 
This may either be in form of huge AGN outbursts which heat 
the surrounding cluster gas, the effect of which lasts for sev- 
eral cooling cycles (such as Hydra-A, MKW3S, A2597 and 
A3581), or in form of short-lived AGN outbur sts which are 



repeat ed after short intervals. A recent study by IShabala et al 



(120081) has shown that the CCRSs in more massive clusters un- 
dergo AGN outbursts more frequently than the AGN in their 
less massive counterparts. Addit ionally, radio source models 
employed bv lShabala et al.l (|2Q08) show that the duration of the 
on-state of an AGN has the same relation with the stellar mass 



as the mass deposition rate has with the stellar mass dBest et al 
l2005h . suggesting the switching on and off of an AGN resulting 
directly from either availabilit y or dep letion of cool cluster gas. 

More recently, IVoit et al.1 (120081) have provided evidence 
that the AGN activity, ICM cooling and star-formation might 
all be linked together through the process of electron ther- 
mal conduction. According to their study, the efficacy of ther- 
mal conduction depends on the size of the temperature in- 
homogeneities relative to the critical length scale associated 
with conduction, Xf. The state of equilibrium between radia- 
tion losses and conduction gain can be equivalently expressed 
in terms of /If and K. Based on above arguments, conduction 
sets an entropy threshold such that only those clusters whose 
central entropy is less than 30 ke V cirT 2 show star formation, in 
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the form of Ha, and enhanced AGN activity. Above this thresh- 
old, conduction is a viable heating mechanism. Further, since 
gas entropy is very closely linked to the central cooling time 
through the relation, f coo i oc K 3/2 /T, such that 0.6 Gyr (1 Gyr) 
corresponds approximately to 30 keV cm~ 2 (43 keV cm -2 ), 
these observations are also in concordance with our results 
on the central temperature drop and cuspiness displayed in 
Figure The central entropy (or central cooling time) thresh- 
old may be an explanation for observing an abrupt central tem- 
perature drop and an increase in cuspiness for clusters with 
cooling times shorter than 1 Gyr. In the cluster regime with 
fcooi < 1 Gyr, AGN heating is the dominant balancing mecha- 
nism to cooling. 

We point out that there is non-negligible scatter in Figures|8] 
and [10] the origin of which could be either extrinsic or intrin- 
sic. In the presence of an AGN-regulated feedback, an intrin- 
sic scatter may imply that the synchrotron luminosity, which 
is only a small fraction of the total AGN output, is not a very 
reliable quantity to use to establish the AGN-ICM interaction. 
It has been noted in previous studies, that the ratio of kinetic 
to radiative AGN power indee d shows a broad range, from a 



few to a several thousands (e.g. lBirzan et al.l2 008: Birz an et al 



2004). Kinetic AGN luminosity may be a more robust measure 
of the total AGN feedback. For this, one requires the radio mor- 
phology of jets and lobes overlaid on X-ray images to help find 
or confirm the X-ray cavities. The census of X-ray cavities is 
highly incomplete since they are of very low contrast, yet they 
are important contributions to the heating budget. An extrin- 
sic scatter in the plots would point at inaccurate measurements 
of observable parameters at both the wavebands, radio (incom- 
plete spectral information) as well as X-ray (imprecise mass 
deposition rates due to spectral resolution power of ACIS on 
Chandra). 

However, strong correlations found in this work between 
the total AGN radio power and various cluster parameters lend 
confidence in synchrotron luminosity as a fairly good measure 
of the cooling activity in clusters. These correlations also pro- 
vide us with motivation to continue our work to acquire low- 
frequency radio measurements for CCRSs which have no reli- 
able data below 500 MHz (constituting 35% of the HIFLUGCS 
sample, see Section l2.1.2b . To achieve this goal, we are await- 
ing proprietary data for all but two of the remaining 35% clus- 
ters with VIA at 325 GHz and with GMRT at dual-frequency 
band 610 MHz/235 MHz. 

5. Conclusions 

We have presented a detailed joint analysis of the brightest 
complete local sample of galaxy clusters, HIFLUGCS, using 
high-resolution X-ray data acquired from Chandra and radio 
data compiled from various sources spanning a wide range of 
frequencies. This study was conducted so as to explore the role 
of AGN in the centers of galaxy clusters accompanied with the 
cooling flow problem. The main results of this study are: 

1 . We find that the integrated radio luminosity (Lr) of a clus- 
ter central radio source is tightly correlated to its 1 .4 GHz 
monochromatic luminosity (see Figure [3) with the excep- 



tion of a few outliers. This correlation is quantified as 
Lr oc Ly}^P l (Equation |2). To estimate the integrated 
radio luminosity of a CCRS, a special effort was made to 
compile low-frequency radio measurements in order to get 
an accurate measure of the total radiative output of the cen- 
trally located AGN. 

2. The best property to diagnose a cool-core cluster with high 
quality data is the central cooling time, f coo i. Based on f coo i, 
there is an increasing probability for the brightest clus- 
ter galaxy (BCG) closest to the X-ray peak emission to 
harbor an AGN with decreasing cooling time (Figure [6). 
The percentage of AGN in three bins ordered in cooling 
time, strong cool core (SCC, f coo i < 1 Gyr), weak cool- 
core (WCC, 1 Gyr < f C00 ] < 7.7 Gyr) and non-cool- 
core (NCC, f coo i > 7.7 Gyr) clusters, is 100%, 67% and 
45% respectively. 

3. There is a trend between the Lr and f coo i such that the for- 
mer increases with decreasing cooling time. This is shown 
in Figure|7] although with a large scatter especially towards 
short f CO oi where the trend appears to break down. 

4. The total radio output of a CCRS scales with the cluster 
size (e.g. X-ray luminosity). This correlation is particularly 
noticeable in SCC clusters (see Figure [8] left panel). The 
best-fit powerlaw for the SCC clusters is L R oc ^ 38±016 
(Equation|6]l. 

5. The total radio output of cool-core clusters (SCC and WCC 
clusters) shows a tight correlation with the classical mass 
deposition rate, M c iasskai (Figure[8| right panel). The corre- 
lation is given by Lr oc ^ c \ a ^f c ° al 25 (Equation |7J. 

6. The radio luminosity of the central radio source shows a 
weak trend with the mass of the supermassive black hole, 
but this trend is seen only for the SCC clusters (FigurelTOl. 
This trend is approximately given by Lr oc Mg^ * 42 . 

7. The NIR bulge luminosity of the BCG (closest to the X- 
ray peak emission) shows a correlation at an unprecedented 
level with the global cluster properties, such as X-ray mass 
shown in the left panel of Figure [TT] and luminosity shown 
in the right panel of Fi gure [TT1 (Equations [TOl and [T2l . 

While we have provided strong evidence of there being an 
abundance of CCRSs with enhanced radio luminosities in clus- 
ters where cooling activity is at its full thrust (SCC clusters), a 
feedback heating-cooling loop may require involvement of ad- 
ditional physical processes, such as conduction as mentioned 
above, to halt the cooling in WCC clusters. 
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Table 1. Radio data details for the central radio sources in the H1FLUGCS sample of galaxy clusters. 
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Redshift 


RA (J2000) 


Dec (J2000) 


v (GHz) 


S v (mJy) 
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Table 1. continued. 



Cluster 


Redshift 


RA (J2000) 


Dec (J2000) 


v (GHz) 


S v (mJy) 
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Table 1. continued. 



Cluster 


Redshift 


RA (J2000) 


Dec (J2000) 


v (GHz) 


S v (mJy) 
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S.I. 
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1080.0+40.0 

2067.0+104.0 
3700.0+130.0 
6600.0+900.0 
(1.04+0.04)E+04 
(2.03±0.25)E+04 


3.990000+0.492000 


3.840000!™ 


-0.95 


6 

6 
6 
6 
6 
6 


A3395s 


0.051000 


06 /! 26'49.74" 
06''26'49.58" 


-54°32'33.6" 
-54°32'34.0" 


4.85 

0.843 
0.408 


1021+54.0 

5610.0+500.0 
7580.0+330.0 


2.050000+0.093300 


1.350000!°;™ 


cs 


6 

6 
6 


A0576" 


0.038900 


07*21' 30.26" 

07 A 2l'30.23" 
07 /, 2l'30.27" 


+55°45'50.6" 

+55°45'41.6" 
+55°45'41.5" 


1.4 


1.6+0.2 


0.000547+0.000068 


0.000582!^ 


-1.0 


1 
4 


A0754 


0.054200 


09 A 09' 16.66" 
09' , 08'32.38" 


-09°4l'20.8" 
-09°37'47.0" 






0.001020 


0.001100 




NO CRS 


A07R0 


O^QOO 
U.UJJ7UU 


oq'm 806 oq" 

\Jy lo \J\j.\)y 

09 /, 18'05.65" 
09 A 18'05.78" 


1 ?°os'4.s o" 

-12°05'43.9" 
-12°05'41.3" 


1 7 

5.0 
1.4 
0.074 
0.0126 


fiQRO 0+1 SO O 

(1.39+0.04)E+04 
(4.51+0.04)E+04 
(57.96+5. 82)E+04 
(390.0±62.4)E+04 


^O 1 0000+0 7SOOOO 


?q ooooo + i-78oooo 

2V - UUUUU -1. 790000 


O QS, 


o 

6 
6 
6 
6 


A 1060 


0.012600 


10 /, 36'42.75" 
10 A 36'42.82" 


-27°3l'42.0" 
-27°3l'42.0" 






0.000052 


0.000053 


- 


NO CRS 


A1367 


0.022000 


11*45' 00.29" 
ll' , 44'02.17" 


+ 19°40'30.2" 
+ 19°56'59.3" 


- 


- 


0.000160 


0.000167 


- 


NO CRS 


MKW4 


0.020000 


12 A 04'27.08" 

12 /, 04'27.05" 
12''04'27.08" 


+01°53'46.1" 

+01°53'45.6" 
+01°53'45.l" 


4.86 
1.4 


0.35+0.1 
2.40+0.5 


0.000213+0.000045 


0-°00815« 


-1.55 


1 
1 



Table 1. continued. 



Cluster 


Redshift 


RA (J2000) 


Dec (J2000) 


v (GHz) 


S v (mJy) 


L v ci.4 GHz> 


L tot 


S.I. 


Ref. 


ZwC11215 


0.075000 


12 A 17'41.7l" 
12''17'41.12" 


+03°39'18.4" 
+03°39'20.8" 






0.004790 


0.005270 




NO CRS 


NGC4636 


0.003700 


12*42'49.9l" 

12''42'49.86" 
12*42'50.0l" 


+02°41'12.6" 

+02°41'16.0" 
+02°41'16.2" 


4.85 

1.4 
0.074 


40.0+3.0 

96.2+1.0 
900.0+160.0 


0.000282±0.000003 


000250+ 0000026 


-0.74 


1 

2 
3 


A3526 


0.011400 


12 A 48'48.85" 

12*48'49.27" 
12''48'49.18" 


-41°18'43.8" 

-41°18'39.9" 
-41°18'42.1" 


5.0 

1.4 

0.843 

0.080 


1320.0+132.0 

3980.0+111.0 
6250.0±40.0 
(3.9±0.6)E+04 


0.H2000+0.003050 


094300+ 0010500 

U.U5"tJUU 001230() 


-0.85 


6 

6 
5 
5 


A 1644 


0.047300 


12*57' 11. 79" 

12*57' 11.57" 
12 A 57' 11.59" 


-17°24'32.3" 

-17°24'34.4" 
-17°24'34.2" 


8.4 

4.85 
1.4 


70.1+0.2 

131.7+0.1 
98.4+3.0 


0.047600±0.00l5l0 


062100+ 0000677 


CS 


1 

1 

2 


A1650 


0.083838 


12*58'41.48" 
12''58'41.49" 


-01°45'42.7" 
-01 o 45'41.0" 






0.000849 


0.000942 


- 


NO CRS 


A1651 


0.084945 


12 /, 59'22.16" 

12 ;, 59'22.51" 
12 /, 59'22.24" 


-04° 11 '49.2" 

-04° 11 '46.0" 
-04° 11 '45.0" 


1.4 

0.324 


7.7+1.1 
22.2+2.0 


0.013200+0.001790 


012100+ 0004200 

—0.003100 


-0.72 


2 
1 


A 1 




iz jy jj. / j 

12*59'35.70" 
12''59'35.47" 


x97°^7^zL Q" 

+z / j / j^+.y 

+27°57'33.8" 
+27°57'36.4" 


2.38 
1.4 
0.408 
0.074 


132.0+6.0 
207.6+1.0 
507.0+72.0 
1500.0+190.0 




U.U2U4UU_ 0003440 


-u. / 1 




6 
2 
6 
6 


NGC5044 


0.008700 


13*15' 23.88" 

13*15'23.96" 
13*15'23.97" 


-16°23'06.8" 

-16°23'07.9" 
-16°23'07.7" 


4.86 
1.4 


28.7+0.2 
36.65+1.0 


0.000595+0.000000 


000701+ 0000019 


-0.20 


1 
2 


A1736 


0.045800 


13 /, 26'51.87" 
13 A 27'28.04" 


-27°10'26.8" 
-27°19'28.8" 






0.000719 


0.000770 




NO CRS 



Table 1. continued. 



Cluster 


Redshift 


RA (J2000) 


Dec (J2000) 


v (GHz) 


S v (mJy) 


L v (1.4 GHz) 


Aot 


S.I. 


Ref. 


A3558 


0.048000 


13 A 27'56.89" 

13''27'56.88" 
\y i 21 56.85 


-31°29'43.2" 

-31°29'43.7" 
-31 29 43.8 


4.86 
1.4 


1.2±0.2 
4.0±0.6 


0.002110+0.000329 


0.002180!°;™ 


-0.97 


1 
2 


A3562" 


0.049000 


13 ;, 33'37.29" 

13 A 33'34.73" 
13 ; '33'34.76" 


-31 o 40'l7.0" 

-31°40'20.1" 
-31°40'20.l" 


1.4 


0.51+0.09 


0.000281+0.000050 


0.000302!°;™ 


-1.0 


12 


A3571 


0.039100 


13 A 47'28.32" 

13 ;, 47'28.38" 
13*47' 28.50" 


-32°51'57.5" 

-32°51'54.0" 
-32°5l'54.0" 


2.38 
1.38 


2.52+0.15 
4.05+0.15 


0.001390+0.000056 


001290+ 0000554 


-0.87 


13 
13 


A 1 7CK 


U.UOZH- / u 


1 ^' ! /t8'S9 SS" 

lj HO JZ.JO 

13 /, 48'52.51" 
13 /, 48'52.43" 


tZO jj 

+26°35'34.8" 
+26°35'33.6" 


A 

■+.OJ 

1.4 
0.408 
0.074 
0.026 


ZD 1 .\J± J^+.KJ 

925.0+27.8 
3150.0+250.0 
(1.09+0.1 l)E+04 
(3.4±0.8)E+04 


U . O^+J WWW ± W . UZ'H-H-U W 


U./88UUU _ 109000 


n Q6 

-\J. VD 


D 

6 
6 
3 
6 


A3581 


0.023000 


14 /, 07'30.19" 

14 ; '07'29.78" 
14 ;, 07'29.82" 


-27°0l'l0.5" 

-27°0l'04.3" 
-27°0l'04.5" 


4.85 

1.4 

0.408 

0.074 


533.0+30.0 

646.2+23.0 
1090.0+40.0 
3310.0+410.0 


0.074500+0.002600 


101000 +000898 ° 

u.iuiuuu _ 008490 


cs 


6 

6 
6 
3 


MKW8 


0.027000 


14''40'43.08" 

14 /, 40'42.87" 
14"40 42.89 


+03°27'57.7" 

+03°27'55.5" 
+03 27 55.2 


4.86 
1.4 


2.09+0. 15 
2.54+0.1 


0.000402+0.000025 


000505+ 0000068 


-0.16 


1 
4 


RXJ1504 


0.215300 


15 A 04'07.52" 

15''04'07.52" 
15 ;, 04'07.52" 


-02°48'16.8" 

-02°48'l6.l" 
-02°48'16.6" 


4.86 
1.4 


46.7+0.1 
67.0+2.0 


0.774000+0.026700 


994000+ 023900 

u.?5"+uuu .164000 


-0.29 


1 

2 


A2029 


0.077280 


15''10'56.06" 

15''10'56.10" 
15''10'55.87" 


+05°44'41.4" 

+05°44'41.6" 
+05°44'39.l" 


4.85 

1.4 
0.408 
0.178 
0.074 


84.0+13.0 

725.0+175.0 
3740.0+170.0 
7000.0+875.0 
(1. 677+0. 172)E+04 


1.060000+0.241000 


1 A ~> AOAAA 

50000 +0 - 29200() 

1.1JUUUU_0.340000 


CS 


6 

6 
6 
6 
6 



Table 1. continued. 



Cluster 


Redshift 


RA (J2000) 


Dec (J2000) 


v (GHz) 


S v (mJy) 


L v (1.4 GHz) 


hot 


S.I. 


Ref. 


A90S9 




1 J ID HO. J 1 

15*16'44.48" 
15 A 16'44.56" 


tu/ yji 17.0 

+07°0l'l8.0" 
+07°01'17.9" 


o.yj 

1.4 

0.408 
0.074 
0.026 


5460.0+100.0 
(2.518±0.109)E+04 
(12.165±1.249)E+04 
(23.2±1.5)E+04 


1 ^0000+0 097800 


1 S70000+ - 110000 
1 ■ 3 /UUUU -0. 108000 






6 
6 
6 
6 


MKW3S 


0.045000 


15*21'51.75" 
15*21'51.87" 


+07°42'28.7" 
+07°42'31.9" 


1.28 

0.61 

0.235 

0.074 


139.0+0.030 

1010.0+0.15 
8436+2.0 
(5.835+0.596)E+04 


0.067800+0.000014 


1 160000+ - 217000 

i.iuuuuu _ 859 o o 


cs 


14 

14 
14 

3 


A2065 


0.072600 


15 A 22'29.32" 

15 A 22'29.16" 
15 22 29.17 


+27°42'22.2" 

+27°42'27.3" 
+27 42 27.5 


4.86 
1.4 


1.2+0.2 
9.2+1.6 


0.012000+0.002220 


026400+ 000992 


-1.50 


1 

3 


A2063 


0.034937 


15*23'05.1l" 

15*23'05.30" 
15*23'05.30" 


+08°36'26.9" 

+08°36'33.0" 
+08°36'33.0" 


4.86 
1.4 


5.2+1.0 
15.7+1.0 


0.004290+0.000282 


004050+ 0002390 


-0.89 


1 
2 


A2142 


0.090900 


15''58'20.65" 

15 /, 58'20.02" 
15"58 20.52 


+27°13'49.2" 

+27° 14' 00.0" 
+27 14 16.5 


1.4 
0.327 


2.5+0.4 
3.2+0.7 


0.004600+0.000708 


005950 +0 004150 

u.uuj^ju _ 002350 


0.18 


1 

11 


A2147 


0.035000 


16 A 02' 16.78" 

16''02' 17.00" 
16 A 02' 16.88" 


+ 15°58'25.6" 

+ 15°58'28.2" 
+ 15°58'28.l" 


4.86 
1.4 


7.9+0.4 
16.1 + 1.1 


0.004370+0.000307 


003900+ 000578 


-0.57 


1 
2 


A2163 


0.203000 


16''15'46.69" 
16''15'48.99" 


-06°09'00.3" 
-06°08'41.5" 






0.017400 


0.021500 




NO CRS 


A2199 


0.030151 


16 A 28'38.32" 

16 /, 28'38.27" 
16''28'38.34" 


+39°33'01.2" 

+39°33'04.9" 
+39°32'60.0" 


5.0 

1.4 

0.750 
0.178 
0.074 
0.010 


487.0+30.0 

3580.0+80.0 
9190.0+120.0 
(5.11±0.26)E+04 
(1 1.24+1. 13)E+04 
(40.8+7. 8)E+04 


0.737000+0.015700 


0.862000~ 


CS 


6 

6 
6 
6 
6 
6 



Table 1. continued. 



Cluster 


Keusniit 


KA (JzUUUj 


Dec (J/UUU) 


V (unz ) 


Jy (mJy) 


L v (i.4 GHz) 


Mot 


C T 


Ker. 


A2204 


0.152158 


16''32'46.94" 

16''32'46.83" 
16 32 46.87 


+05°34'31.3" 

+05°34'31.8" 
+05 34 34.9 


1.4 

0.365 
0.074 


69.3+2.5 

276.0+27.0 
1640.0+190.0 


0.427000+0.014500 


576000+ - 106000 
u.j/uuuu_ 184000 


-1.08 


2 

4 
3 


A2244 


0.096800 


17''02'42.68" 

17''02'42.47" 
17''02'42.3l" 


+34°03'39.3" 

+34°03'36.3" 
+34°03'38.2" 


1.46 
0.327 


3.1+0.4 
5.8+1.0 


0.006770+0.000823 


007120+ 0003210 
u.uu/ iz.u_ 002150 


-0.42 


2 
1 


A2256 


0.058100 


17''03' 14.26" 
17*03'35.9l" 


+78°38'59.9" 
+78°37'43.5" 


- 


- 


0.00 ll 80 


0.001280 


- 


NO CRS 


A2255 


0.080600 


17''12'34.15" 
17*12'28.75" 


+64°04'11.5" 
+64°03'38.5" 


- 


- 


0.002340 


0.002590 


- 


NO CRS 


A3667 


0.055600 


20 /l 12'42.66" 
20''12'27.26" 


-56°50'48.6" 
-56°49'36.3" 


- 


- 


0.00 1 290 


0.001400 


- 


NO CRS 


S1101 


0.058000 


23''13'58.40" 

23 /, 13'58.63" 
23 ;, 13'58.63" 


-42°43'31.0" 

-42°43'39.3" 
-42°43'39.9" 


4.85 

0.843 
0.408 


49.0+10.0 

473.56+10.0 
1390.0+60.0 


0.192000+0.020300 


595000 +015900 ° 


-1.39 


6 

5 
6 


A2589 


0.041400 


23 /, 23'57.40" 
23 A 23'57.4l" 


+ 16°46'37.9" 
+ 16°46'37.9" 






0.000583 


0.000622 


- 


NO CRS 


A2597 


0.085200 


23 /, 25' 19.93" 

23 ;, 25' 19.73" 
23*25' 19.82" 


-12°07'27.5" 

-12°07'27.5" 
-12°07'28.6" 


4.85 

1.4 

0.408 

0.074 


407.0+24.0 

1875.0+56.0 
7200.0+400.0 
(2.87+0.29)E+04 


3.320000+0.093300 


3 060000 +a291000 

j.uuuuuu _ 0580000 


CS 


6 

6 
6 
6 


A2634 


0.031385 


23''38'29.25" 

23 /, 38'29.38" 
23 /, 38'29.50" 


+27°01'54.2" 

+27°01'52.6" 
+27°01'53.8" 


5.0 

1.4 
0.074 


2800.0+80.0 

7900.0+400.0 
(6.95±0.37)E+04 


l. 730000+0.086500 


1 500000+ 0111000 

i.juuuuu Q 106000 


-0.86 


6 

6 
3 



Table 1. continued. 



Cluster 


Redshift 


RA (J2000) 


Dec (J2000) 


v (GHz) 


S v (mJy) 


L v (1.4 GHz> 


L tot 


S.I. 


Ref. 


A2657 


0.040200 


23 ;, 44'57.48" 
23''44'57.42" 


+09°ll'31.0" 
+09°ll'34.9" 






0.000549 


0.000585 




NO CRS 


A4038 


0.030000 


23 /, 47'43.18" 

23''47'45.04" 
23''47'44.95" 


-28°08'31.2" 

-28°08'26.5" 
-28°08'24.5" 


4.86 

1.4 
0.074 


9.4±0.1 

25.0+0.030 
(1. 375+0. 142)E+04 


0.005190+0.000008 


190000+ 0027900 

u.ijuuuUj 164000 


cs 


1 

7 
3 


A4059 


0.047500 


23''57'00.93" 

23 /, 57'00.68" 
23*57' 0.67" 


-34°45'33.3" 

-34°45'33.l" 
-34°45'31.7" 


4.85 

1.4 

0.408 

0.074 


117.0+13 

1286.0+44 
8700.0+350.0 
(5.762+0.589)E+04 


0.682000+0.022500 


1 150000 +a 148000 


cs 


6 

6 
6 
3 



Columns - (1) Cluster name, (2) Redshifts compiled from NASA/IPAC Extragalactic Database (NED), (3) & (4) Right Ascension and Declination: The 
first row corresponds to the X-ray peak position {Chandra), the second row corresponds to the optical position (2MASS) and the third row corresponds to 
the radio position (at the highest available resolution), (5) Radio frequency in GHz, (6) Radio flux density in mJy, (7) Spectral radio luminosity at 1 .4 GHz 
in units of 10 32 erg s _1 Hz -1 , (8) Total radio luminosity between 10 MHz and 15 GHz in units of 10 42 erg s , (9) SI: Spectral index of a radio source. 
'CS' denotes a complicated spectrum with breaks and/or turnovers, (10) References for the flux density. 

References - (1) This work - based on data from the VLA (Very Large Array) archive (2)NVSS- The NRAO VLA Sky Survey 
at 1.4 GHz [Condon et al. 1998] (3) VLSS - The VLA Low-Frequency Sky Survey at 74 MHz [Cohen et al. 2007] (4) The VLA 
FIRST (Faint Images of the Radio Sky at Twenty-cm) Survey at 1.4 GHz (5) SUMSS - The Sydney University Molonglo Sky Survey at 
843 MHz [Bock et al. 1999] (6) NED - NASA/IPAC Extragalactic Database (7) Slee et al. 2001 (8) Birzan et al. 2004 (9) Blanton et al. 2004 
(10) Patnaik& Singh 1988 (11) Markovic et al. 2004 (12) Miller et al. 2005 (13) Venturi et al. 2002 (14) Giacintucci et al. 2007 
'NO CRS': Clusters with no central radio source within 50 kpc of the X-ray peak. We calculated L x 4 qh z and Lr based on 3-cr upper-limit, where <x 
is the background rms in the map, and assuming S oc y~ 10 . Except two clusters for which we used data from the VLA archive, the upper-limits are based 
on either NVSS or MOST observations. 

a A576, A3 158 and A3562 have measurements only at a single frequency; hence we assume S oc y~ 10 . 

b A3376 is an exceptional case where the 'brightest' cluster galaxy (km ext = 10.38) is located ~ 1 Mpc from the X-ray peak. However, for the purpose of 
our study we assumed the powerful wide-angled-tail radio galaxy (km ext = 11.11) very close to the X-ray peak (< 15 kpc) as the "BCG". The origin of 
radio activity is disputable; it could either be related to ICM cooling, the small displacement from the X-ray peak having being caused by a merger, or it 
could have been triggered by the merger itself. 



